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ABSTRACT 


This  thesis  considers  the  problem  of  generating  optimal  entry  trajectories  for  a 
reusable  launch  vehicle  following  control  surface  failures.  The  thesis  builds  upon  the 
work  of  Dr.  David  Doman,  Dr.  Michael  Oppenheimer  and  Dr.  Michael  Bolender  of  the 
Air  Vehicles  Directorate,  Air  Force  Research  Lab  Dayton  Ohio.  The  primary  focus  of 
this  work  is  to  demonstrate  the  feasibility  of  inner  loop  reconfiguration  and  outer  loop 
trajectory  retargeting  and  replanning  for  the  X-33  reusable  launch  vehicle  (RLV) 
following  the  imposition  of  a  control  surface  failure.  The  trajectory  generation  model 
employs  path  constraints  generated  by  an  AFRL  trim  deficiency  algorithm  coupled  with 
an  inner  loop  control  allocator  and  aerodynamic  database  that  captures  the  full  6-DOF 
vehicle  aerodynamic  effects  while  utilizing  an  outer  loop  3-DOF  model.  The  resulting 
optimal  trajectory  does  not  violate  the  trim  deficiency  constraints  and  provides  additional 
margins  for  trajectories  flown  during  failure  conditions.  The  footprints  generated  by  the 
thesis  show  that  contemporary  footprint  analysis  for  vehicles  experiencing  control  surface 
failures  are  overly  optimistic  when  compared  to  those  footprints  that  consider  vehicle 
aerodynamic  stability  and  realistic  landable  attitudes  at  the  threshold  of  the  landing 
runway.  The  results  of  the  thesis  also  show  the  performance  reductions  resulting  from 
decoupling  the  inner  and  outer  loop  and  that  trajectories  can  be  generated  to  the  landing 
runway  without  using  a  region  of  terminal  area  energy  management. 
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I.  INTRODUCTION 


A.  HISTORICAL  BACKGROUND 

As  early  as  1931,  Dr.  Eugen  Sanger  eoneeived  of  an  aerovehiele  that  eould 
aehieve  low-earth  orbit  (LEO)  and  return  to  earth  like  a  glider.  The  Sanger  Silverbird 
was  advoeated  as  a  hypersonie  vehiele  that  eould  deliver  an  8300  pound  payload  into 
LEO  and  return  to  earth  in  an  unpowered  glide.  Alternatively,  the  Silverbird  (or  its 
military  inearnation,  the  Amerika  Bomber)  eould  deploy  a  1600-pound  payload  into  a 
suborbital  trajeetory  from  a  launeh  point  halfway  around  the  globe  and  return  to  its 
origination  by  skip-gliding  or  dynamie  soaring  [ref  25]. 

The  problem  of  extreme  reentry  heating  was  severely  underestimated  by  Sanger 
and  eonelusions  from  NACA  Ames  reentry  studies  in  the  1950’s  led  to  the  abandonment 
of  skip  gliding  or  dynamie  soaring  beeause  there  was  no  eontemporary  teehnologieal 
solution  to  mitigate  the  severe  heating  during  reentry.  Later  manned  missions  sueh  as 
Apollo  overeame  the  reentry  heating  problem  by  employing  a  vehiele  that  used  an 
ablative  heat  shield  to  proteet  the  payload  and  followed  a  steep  ballistie  reentry  trajeetory. 
Unfortunately  this  type  of  teehnology,  while  effeetive,  did  not  allow  for  reusable 
vehieles. 

Interest  in  reusable  spaeeeraft  resulted  in  the  development  of  the  spaee  shuttle 
whieh  makes  use  of  advaneed  materials  and  eontrols  teehnologies  to  fly  trajeetories 
whieh  keep  the  shuttle  within  a  very  narrow  reentry  eorridor  to  limit  the  heating  and 
dynamie  pressure  effeets  of  the  entry  trajeetory  thereby  ensuring  a  sueeessful  landing 
[ref.  22].  These  reentry  trajeetories  are  eomputed  offline  for  a  variety  of  initial 
eonditions  and  then  stored  in  the  shuttle’s  eomputer  where  they  ean  then  be  aeeessed  by 
the  guidanee  eomputer.  The  generation  of  these  trajeetories  represents  a  signifieant  effort 
and  is  supported  by  a  large  body  of  work  from  eontemporary  optimal  eontrol  researeh. 

Reusable  launeh  vehieles  are  being  pursued  as  a  low-eost  alternative  to 
expendable  launeh  vehieles  and  the  shuttle.  The  eoneept  of  a  reusable  launeh  vehiele  that 
ean  aehieve  LEO  using  a  single  stage  is  seen  as  a  teehnology  that  ean  signifleantly  reduee 
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the  cost  of  access  to  space  as  well  as  improving  the  frequency  of  access  to  space.  The 
employment  of  autonomous,  reusable  launch  vehicles  requires  additional  control 
robustness  to  replace  the  human’s  place  in  the  loop.  The  guidance  and  control  of  these 
autonomous  reusable  launch  vehicles  mandate  new  guidance  strategies  that  are  able  to 
identify  and  adapt  to  vehicle  failures  during  the  trajectory  and  still  return  to  earth  safely. 

An  important  class  of  failures  that  concern  reentry  trajectories  are  those  control 
surface  failures  which  affect  the  vehicle’s  aerodynamic  characteristics  and  controllability. 
The  nominal  reentry  trajectories  are  designed  assuming  the  vehicle’s  original 
aerodynamic  characteristics  and  can  no  longer  be  assumed  valid  following  a  control 
surface  failure.  The  challenge  is  to  reconfigure  the  inner  loop  control  algorithm, 
determine  the  new  aerodynamic  characteristics  of  the  failed  vehicle  and  to  generate  a  new 
trajectory  using  the  new  aerodynamic  data  to  direct  the  vehicle  to  a  landing  point  using  a 
modified  performance  criterion. 

This  thesis  demonstrates  the  application  of  a  psuedospectral  method  to  determine 
the  optimal  entry  trajectory  for  the  X-33  for  nominal  conditions  as  well  as  for  off- 
nominal  conditions  caused  by  vehicle  control  surface  failures.  Actual  X-33  aerodynamic 
data  is  used  in  a  reduced-order  dynamical  model  to  generate  the  trajectories.  Feasibilty, 
optimality  and  various  performance  metrics  of  the  psuedospectral  approach  are 
investigated  and  quantified. 


B,  THE  X-33  PROGRAM 

The  X-33  was  a  half-scale  technology  demonstrator  of  a  potential  future  Reusable 
Launch  Vehicle  (RLV)  [ref  24].  The  prime  contractor  for  the  1.2  billion  dollar  program 
was  Lockheed  Martin  who  called  the  full-scale  vehicle  VentureStar.  The  X-33  was  to 
provide  data  that  would  lead  to  the  full  development  of  a  single  stage  to  orbit  (SSTO) 
vehicle  that  would  dramatically  lower  the  cost  of  access  to  space  through  SSTO  and 
minimal  time  between  missions.  The  X-33  represented  cutting  edge  technology  in  all 
design  aspects,  from  its  linear  aerospike  engines  to  it  composite  fuel  tanks.  The  vehicle 
can  reach  altitudes  of  50  miles  and  hypersonic  speeds.  The  design  turn  around  time  for 

this  vehicle  was  7  days  with  an  eventual  turn  around  goal  of  two  days.  The  program  was 
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expected  to  begin  full  flight  demonstration  in  1999  but  was  discontinued  in  2001  after 
technology  maturation  difficulties.  The  X-33  remains  a  desirable  vehicle  to  study 
because  of  the  large  amount  of  performance  data  available  and  because  of  its  multiple 
redundant  control  surfaces  allows  for  the  investigation  of  a  considerable  variety  of 
recoverable  control  surface  failures  and  their  impacts  upon  entry  trajectories. 

1,  X-33  Mission  Description 

The  X-33  was  designed  to  demonstrate  the  ability  to  carry  a  small  payload  to  low 
earth  orbit,  return  safely  to  earth,  and  undergo  minimal  preparations  in  order  to  be  ready 
for  another  mission.  The  vehicle  was  to  be  launched  from  Edwards  Air  Force  Base  and  is 
potentially  recoverable  at  a  number  of  suitable  fields.  The  proposed  flight  test  program 
utilized  several  ascent  trajectory  profiles  terminating  at  Malstrom  Air  Force  Base, 
Montana  or  Michael  Army  Air  Field,  Fltah. 

The  typical  X-33  mission  begins  with  a  powered  ascent  that  achieves  a  given 
Mach  and  altitude  at  main  engine  cutoff  (MECO).  The  vehicle  performs  some  activity  on 
orbit,  transitions  and  enters  the  atmosphere  passing  through  several  intermediate  flight 
stages  ultimately  ending  at  landing  and  rollout.  Figure  1  details  the  various  X-33 
mission  phases  as  well  as  the  effectors  that  are  used  to  control  the  vehicle.  While  Figure 
1  shows  a  reaction  control  system  (RCS)  available  during  the  entry  phase,  this  thesis 
assumes  that  no  RCS  present  during  the  entry  portion  of  the  mission.  The  weight  on 
wheels  (WOW)  and  rollout  mission  phase  are  not  considered  in  the  trajectory  generation. 
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Figure  1  X-33  Flight  Phases  (from  Ref.  5) 


Figure  2  depicts  typical  X-33  mission  profiles  with  recoveries  at  both  fields. 
Notice  that  the  ascent  profiles  and  the  entry  initial  conditions  differ  significantly  for  each 
mission. 
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Figure  2  Representative  X-33  Mission  Profiles  (from  Ref.  5) 

A  partial  list  of  mission  profiles  is  shown  in  Table  1.  The  Maeh  number  and 
maximum  altitude  are  used  as  initial  conditions  for  the  trajectory  optimization  routines. 


Flight 

No. 

Trajectory 

Designation 

Landing 

Site 

Dist 

nm 

Mach 

No. 

Max  Alt. 
k  Feet 

Time 

MECO 

asmg 

Q 

Q 

Alpha 

Fleat 

Rate 

Test  Objectives 

1 

9a-2 

Benign 

Michael 

391 

8.84 

186 

239 

4.11 

196 

4864 

17.27 

Aero  plume  interactions.  Structural 
interachons,  Thrust  Vector  Control 
Verification  Installed  Isp 

2 

9d-l 

Transition 

Michael 

391 

10.66 

186 

197 

4.11 

404 

6286 

24.14 

Shock-Shock  interactions 

Boundary  layer  transition 

3 

7b-1 

Max  Heatinq 

Michael 

391 

8.97 

189 

221 

4.11 

239 

6576 

27.29 

TPS  panel  thermal  properties 

Aircraft-Like  Operational  Efficiency 

4 

9c-2 

Real  Gas 

Michael 

391 

1021 

189 

204 

4.11 

399 

6006 

22.12 

Real  Gas  Effects  Verification 

Aircraft-Like  Operational  Efficiency 
Verification  48  hour  turn-around 

5* 

7b-l 

Max  Heating 

Michael 

391 

8.97 

189 

221 

4.11 

239 

6576 

27.29 

TPS  panel  thermal  propertes  /  support 
structure 

6 

7d 

Real  Gas 

Malmstrom 

826 

12.64 

302 

208 

4.03 

342 

7080 

37.05 

Real  Gas  Effects  Verification 
-  Surface  catalysis  with  atomic  oxygen 

7* 

7a 

Real  Gas 

Malmstrom 

826 

12.64 

302 

208 

4.03 

342 

7080 

37.05 

-  Effects  on  pitch  moment 

8-15 

As  Required 

Michael  / 
Malmstrom 

391/ 

826 

Repeat  Wghts  as  required  to  aquire 
additional  data 

Table  1 .  Flight  Test  Trajectory  Matrix  (from  Ref.  5) 


2,  X-33  Vehicle  Description 

The  X-33  is  a  lifting  body  with  four  sets  of  aerosurfaces:  rudders,  body  flaps,  and 
inboard  and  outboard  elevons.  The  vehicle  also  has  a  reaction  control  system  (RCS)  for 
use  in  unpowered  flight  at  low  dynamic  pressure.  For  powered  flight,  the  vehicle  uses 
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aerospike  engines  with  thrust  veetoring  for  piteh  and  roll  eontrol  and  differential  thrust 
for  yaw  eontrol  [ref  28].  Figure  3  depicts  the  general  X-33  configuration  and  identifies 
the  vehicle’s  control  surfaces. 


Figure  3  X-33  Vehicle  Depiction  (from  Ref  1) 

The  vehicle  gross  liftoff  weight  is  280,000  pounds.  The  vehicle  weight  after 
MECO  is  approximately  78,000  pounds. 

The  control  surfaces  deflection  vector,  also  known  as  the  effector  vector,  is 
identified  by  the  Greek  letter  5  whose  elements  are  identified  as 

J  =  [Relevon,^  Relevon^^^  Lelevon.^  Lelevon^^^  Rudder^j  Rudder^^]^  (1-1) 

The  X-33  control  surfaces  have  maximum  and  minimum  deflections  in  degrees  as 
shown  by  the  vectors  below: 

J^,„=[-30  -30  -30  -30  -15  -15  -60  -30f  (1.2) 

^max=[30  30  30  30  26  26  30  60f  (1.3) 
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Locked  control  surface  failures  are  modeled  by  replacing  elements  of  (1.2)  and 
(1.3)  that  correspond  to  the  failed  surfaces,  with  the  positions  of  the  failed  surfaces. 

C.  IMPLICATIONS  OF  CONTROL  SURFACE  FAILURES 

Many  published  works  deal  with  the  subject  of  reconfigurable  control  following  a 
control  surface  failure  of  the  vehicle  at  some  point  in  its  mission.  The  significance  of 
reconfigurable  control  is  that  it  allows  the  inner  loop  to  compensate  for  the  changed 
aerodynamics  of  the  vehicle  and  maintain  controllability.  For  many  cases,  reconfigurable 
inner  loop  control  is  insufficient  to  continue  safe  flight  and  a  new  trajectory  based  on 
knowledge  of  the  type  of  failure  and  the  resulting  changes  in  the  constraints  on  the  state, 
cost,  control  and  path  arguments  to  the  optimal  trajectory  problem  must  be  computed. 
The  replanned  trajectory  has  the  ability  to  continue  the  mission  from  the  point  of  failure 
and  on  to  landing  thus  saving  the  vehicle  and  crew  which  would  have  otherwise  been  lost 
as  a  result  of  the  failure.  The  availability  of  an  online  algorithm  such  as  the  AFRL  trim 
deficiency  algorithm  provides  the  outer  loop  with  the  feasible  range  of  Mach  and  angle  of 
attack  combinations  throughout  the  trajectory  as  well  as  the  full  6-DOF  effects  of  the  X- 
33  aerodynamic  model  for  use  in  coefficients  of  lift  and  drag  calculations  for  the 
remainder  of  the  trajectory. 

The  thesis  model  uses  a  pseudospectral  method  to  determine  the  optimal  entry 
trajectory  subject  to  appropriate  constraints  such  as  normal  load  and  trim  deficiency.  The 
approach  assumes  online  identification  of  the  failure.  Following  the  failure  identification 
and  characterization,  the  new  effector  vectors  containing  the  upper  and  lower  control 
surface  bounds  are  passed  to  the  trim  deficiency  algorithm  which  determines  the  new 
feasible  ranges  of  Mach  number  and  angle  of  attack  as  well  as  optimized  effector 
displacement  vectors  for  given  Mach-dependent  angle  of  attack  requests.  These  values 
are  passed  to  the  trajectory  algorithm  as  a  path  constraint  to  be  used  while  determining 
the  optimal  trajectory.  An  additional  function  of  the  trim  deficiency  algorithm  queries 
the  vehicle  aerodynamic  database  to  provide  new  vehicle  coefficients  of  lift  and  drag  in 
Mach-alpha  space.  The  outer  loop  dynamical  model  employs  reduced-order  equations  of 
motion  that  use  interpolated  table  lookups  to  determine  the  appropriate  lift  and  drag 
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coefficients  derived  from  the  vehicle  state  (based  on  angle  of  attaek  and  Maeh  number). 
Onee  a  eontrol  surfaee  failure  is  detected  and  identified,  the  aerodynamic  table  is  changed 
to  reflect  the  new  vehiele  aerodynamies  from  the  perspeetive  of  the  reduced  order 
dynamical  model. 
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II.  APPROACH 


A.  GENERAL  APPROACH  DISCUSSION 

The  model  used  in  this  thesis  uses  a  reduced-order  dynamical  model  that  queries 
the  vehicle  full  order  aerodynamic  model  in  order  to  capture  the  full  6-DOF  model 
characteristics.  A  flat  stationary  earth  in  Cartesian  coordinates  is  used  in  order  to  better 
understand  and  visualize  the  interactions  between  the  different  components  of  the 
optimization  routine. 

B,  REDUCED-ORDER  DYNAMICAL  MODEL 

1.  Coordinate  Systems 

The  simulation  uses  a  three-dimensional  Cartesian  local  horizon  coordinate 
system  that  chooses  as  its  origin  the  vehicle's  position  at  the  beginning  of  its  unpowered 
descent  or  at  the  point  in  the  descent  at  which  a  failure  is  identified.  The  positional 
elements  are  expressed  as  pure  states  in  the  easily  recognized  forms  of  x  (downrange),  y 
(crossrange),  and  z  (altitude).  The  coordinate  sign  convention  is  the  usual  xyz  sign 
convention  where  altitude  decreasing  moving  "down",  and  down  range  increases  to  the 
"right"  when  viewing  the  xz  plane  from  the  side.  The  y  axis  is  the  cross  product  of  the  z 
and  X  axes.  Full  treatment  of  dynamical  formulations  in  other  coordinate  systems  can  be 
found  in  Rea  [ref.  7].  Figure  4  shows  the  coordinate  system  used  in  the  reentry 
simulation. 

where  /  =  flight  path  angle,  the  angle  made  by  the  velocity  vector  and  the  x-y 

plane. 

P  =  azimuth  angle,  the  angle  made  by  the  projection  of  the  velocity 
vector  onto  the  Newtonian  x-z  plane  and  the  Newtonian  x-axis.  (This  is  used 
interchangeably  with  heading  angle  since  we  are  dealing  with  a  fiat  earth) 

(j)  =  bank  angle,  the  angle  made  by  the  lift  vector  from  vertical. 

V  =  velocity  magnitude 

a  =  vehicle  angle  of  attack 
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The  body  axes  are  depicted  in  blue;  the  Newtonian  Frame  is  in  black.  The 
velocity  vector  and  its  projections  are  red.  The  lift  and  drag  vectors  are  colored  orange. 
The  vehicle  is  assumed  to  be  in  coordinated  flight  (trimmed  in  yaw,  no  sideslip). 

z 


2,  Equations  of  Motion 

The  equations  of  motion  are  a  simplified  version  of  the  traditional  spherical  earth 
equations  of  motion.  There  are  six  equations  representing  the  kinematics  and  the  3DOF 
dynamics  of  the  problem  as  adapted  from  Wiesel  [ref.  6].  These  equations  of  motion  are 
shown  below; 

a.  Kinematics: 


x  =  V  cos  P  cosy 

(1.4) 

y  =  V  sin  cosy 

(1.5) 

z  =  F  sin  y 

(1.6) 

Dynamics: 

D 

(1.7) 

V  = - g  sm  y 

m 
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(1.8) 


f  = 
0  = 


Lcos^ 

mV 

Lsmcj) 

mV 


gCOSf 

V 


(1.9) 


The  expressions  for  the  lift  and  drag  forces  come  directly  from 
fundamental  aerodynamics,  repeated  here  for  completeness. 

lV-c,a,^pv-  (1.10) 

D  =  \CoA^pV^  (1.11) 

In  a  full  6-DOF  model  the  coefficients  of  lift  and  drag, 
C^[a,M,(7,S),Cjj[a,M,a,S)arQ  functions  of  several  variables:  a  is  the  vehicle’s 

angle  of  attack,  M  is  the  Mach  number,  <j  is  the  sideslip  angle  and  5  is  the  control 

surface  deflection  vector  previously  discussed.  The  reduced  order  model  employed  by 
this  thesis  assumes  that  the  coefficients  of  lift  and  drag  are  functions  of  the  state  variables 
only. 

C^,C,=f{M,a)  (1.12) 

Additionally  Mach  number  and  atmospheric  density  are  functions  of 

altitude: 


M=M{z)  (1.13) 

p  =  p{z)  (1.14) 


This  simplified  assumption  is  valid  because  the  model  performs  a  table 
lookup  of  the  lift  and  drag  coefficients  that  were  obtained  from  a  full  order  model  that 
uses  the  vehicle  flight  condition  and  the  optimized  effector  displacement  vector  as 
arguments  to  interpolate  an  aerodynamic  data  table  for  the  vehicle.  The  resulting 
coefficients  of  lift  and  drag  from  wing,  body  and  trim  effects  are  summed  and 
incorporated  into  the  table.  The  effector  displacement  is  an  implicit  argument  in  the 
algorithm’s  coefficients  of  lift  and  drag  determination.  The  use  of  the  optimized  effector 

vector  effectively  decouples  the  outer  loop  from  the  inner  loop. 
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3,  The  State  Vector 

The  state  veetor  is  ehosen  to  fully  represent  the  state  of  the  system  in  question  at 
any  instant.  Using  the  flat  earth  eoordinate  system,  the  position  of  the  vehiele  ean  be 
represented  by 


X 

T 

z 

where 

X  :=  downrange  distanee  in  feet 

y  :=  erossrange  distanee  in  feet 

z  :=  altitude  in  feet 


(1.15) 


The  veloeity  of  the  vehiele  is  represented  in  spherieal  eoordinates: 

V  :=  the  magnitude  of  the  veloeity  in  ft/see 

Y  :=  the  flight  path  angle  (FPA)  in  radians 

P  :=  the  azimuth  angle,  in  radians 


These  three  state  variables  when  eoupled  with  the  positional  information  eomprise 
the  state  veetor: 


X 

T 


X  = 


z 

V 

r 


P 


(1.16) 


4,  The  Control  Vector 

The  eontrol  veetor  for  this  problem  is  the  veetor  whose  elements  are  the 
eontrollable  parameters  of  the  system.  For  a  vehiele  trimmed  in  yaw,  the  eoeffieients  of 
lift  and  drag  and  therefore  x  and  z  motion,  are  determined  by  the  vehiele's  angle  of  attaek. 
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The  vehicle's  bank  angle  when  combined  with  angle  of  attack  determines  the  vehicle's 
cross  range  motion.  Since  reentry  is  accomplished  dead-stick  there  is  no  thrust  control 
variable.  A  suitable  candidate  for  the  control  vector  is  therefore 


(1.17) 


a  :=  the  vehicle  angle  of  attack 

<j  :=  the  vehicle  bank  angle 

This  type  of  control  is  also  called  inertialess  control  [ref  20] 


C.  4-DOF  MODEL 

The  inertialess  controller  can  be  adapted  to  more  accurately  model  the  physical 
system  by  recognizing  that  physical  actuators  are  rate  limited.  Inertia  can  be  added  to  the 
system  by  imposing  realistic  rate  constraints  on  the  actuators.  This  is  accomplished  by 
incorporating  the  elements  of  the  control  vector  in  Equation  (1.17)  into  the  state  vector  in 
Equation(1.16).  The  new  state  vector  becomes  an  8-element  vector. 

X 

T 
z 
V 

(1.18) 

r 

P 

a 

<j 

The  new  control  vector  components  are  now  the  rates  of  the  previous  control 
vector’s  components  in  Equation(1.17). 


(1.19) 
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This  alternative  formulation  is  provides  a  more  realistie  model  by  preventing 
unlimited  eontrol  surfaee  defleetion  rates.  Note  that  while  the  physieal  eontrols  remain 
a  and  cr ,  the  algorithm  uses  the  rates  of  these  variables  in  its  optimization.  This 
dynamical  model  is  called  a  “4-DOF”  model  because  it  provides  additional  fidelity 
beyond  the  3-DOF  model.  The  full  4-DOF  trajectory  optimization  algorithm  is  shown  in 
Figure  6 


It 


li 


Error 


Optimization  Routine 


Figure  5  Trajectory  Optimization  Feasibility  Error  Estimation 

The  magnitude  of  the  error  shown  in  Eigure  5  is  used  to  choose  the  size  of  the  nodes 
computed  in  the  psuedospectral  method.  The  number  of  nodes  is  increased  until  the  error 
between  the  propagated  states  and  optimal  states  is  acceptable.  Accuracy  comes  at  a 
cost;  the  larger  number  of  nodes  significantly  increases  algorithm  execution  time. 


D,  ATMOSPHERIC  MODEL 

An  AERE  atmospheric  model  which  generates  speed  of  sound  and  air  density  for 
a  given  an  altitude  is  used.  The  outputs  of  the  atmospheric  model  are  combined  with  the 
vehicle  velocity  to  determine  the  vehicle’s  Mach  number  to  be  used  in  the  aerodynamic 
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table  lookups  as  well  as  the  control  allocator  and  trajectory  path  constraints.  Figure  6 
and  Figure  7  depict  the  model’s  variance  of  the  speed  of  sound  and  atmospheric  density 
as  a  function  of  altitude.  The  atmospheric  model  is  a  combination  of  exponential 
functions  separated  by  discrete  altitude  breakpoints. 
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Figure  6  Plot  of  Density  vs.  Altitude  (from  AFRL  atmosphere  model) 
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Figure  7  Plot  of  Speed  of  Sound  vs.  Altitude  (from  AFRL  atmosphere  model) 

E.  GRAVITATIONAL  MODEL 

The  gravitational  model  used  is  a  low  order  model  that  uses  a  constant 
gravitational  acceleration  which  is  assumed  to  always  act  in  the  negative  z  direction.  The 
decision  to  use  a  simplistic  model  is  based  primarily  on  the  desire  to  remove  additional 
complications  and  observe  the  simulation  dynamics  in  as  pure  a  manner  as  possible. 


F.  X-33  AERODYNAMIC  MODEL 

This  thesis  considers  the  atmospheric  entry  of  the  X-33  reusable  launch  vehicle 
following  a  mission  in  low-earth  orbit.  The  computed  trajectory  is  an  unpowered  descent 
to  a  known  landing  field.  The  trajectory  exposes  the  vehicle  to  full  ranges  of  Mach  and 
angle  of  attack.  To  accurately  model  the  aerodynamic  performance  of  the  X-33,  the 
trajectory  code  utilizes  a  dynamic  linked  library  (dll)  developed  by  Dr.  David  Doman  of 
the  Air  Force  Research  Laboratory  (AFRL)  consisting  of  compiled  C++  code  that  queries 
X-33  aerodynamic  data  derived  from  the  full  MAVERIC  simulation.  The  file, 
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x33aero.dll,  determines  the  full  order  aerodynamie  eharaeteristies  ineluding  wing-body 
and  trim  effeets  from  an  input  veetor  eonsisting  of  Maeh  number,  angle  of  attaek  and 
eontrol  effeetor  position.  Sinee  the  position  of  the  eontrol  effeetors  is  arbitrary  for  any 
given  desired  angle  of  attaek,  two-dimensional  tables  of  total  lift  and  drag  eoeffieients  in 
Maeh-alpha  spaee  are  generated  using  another  AFRL  eode  that  uses  a  pieeewise  linear 
programming  (PLP)  algorithm  to  determine  the  optimal  eontrol  alloeation  of  the  vehiele's 
eontrol  surfaees  (i.e.,  eontrol  indueed  moments  balanee  the  moments  generated  by  the 
wing-body  with  minimum  eontrol  defleetion  in  a  1-norm  sense)  to  aehieve  a  desired 
angle  of  attaek.  The  optimal  eontrol  alloeation  is  not  direetly  used  by  the  outer  loop 
trajeetory  generator,  but  the  resulting  eoeffieients  of  lift  and  drag  derived  from  the 
optimized  eontrol  effeetors  for  a  given  eommanded  angle  of  attaek  and  Maeh  are  used  in 
the  outer  loop. 

1.  Mach  Breakpoints 

The  X-33  aerodynamie  data  eovers  a  range  of  angles  of  attaek  in  inerements  of 
one  degree  from  -10  to  50  degrees  over  a  diserete  Maeh  range  given  by  [0  0.3  0.6  0.8  0.9 
0.95  1.05  1.15  1.46  2.01  2.75  3.49  4.75  6  10  20].  This  diserete  Maeh  interval  is  used  in 
the  eompiled  C++  eode  refereneed  by  the  DIDO  dynamies  funetion.  The  aerodynamie 
database  output  is  in  the  body  frame,  so  an  appropriate  transformation  to  a  Newtonian 
frame  is  neeessary  in  order  to  analyze  the  trajeetory.  The  normal  foree  path  eonstraints 
ean  be  analyzed  in  the  body  frame. 

2.  AFRL  PLP  Control  Allocation  Model 

The  AFRL  eontrol  alloeation  algorithm  is  used  to  trim  the  vehiele  as  well  as 
eomputed  the  lift  and  drag  estimates  for  the  lift  in  any  eonfiguration.  The  algorithm 
assumes  that  the  vehiele  is  in  symmetrie  flight  (sideslip  angle  is  zero)  and  the  resulting 
lateral  wing-body  foree  and  moment  eoeffieients  are  zero.  By  assuming  a  zero  magnitude 
body-axis  angular  veloeity,  the  wing-body  pitehing  moment  of  the  vehiele  beeomes  a 
funetion  of  Maeh  number  and  angle  of  attaek.  This  allows  for  the  neeessary  eonditions  to 
rotationally  trim  the  vehiele  that  the  moments  resulting  from  the  eontrol  surfaee 
defieetions  must  be  equal  and  opposite  that  of  those  produeed  by  the  wing-body  [ref.  22] 
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C„yM,a,d) 

= 

CyM,a,d)\ 

0 

Where  l,m,n  denote  the  roll,  pitch  and  yaw  axis,  is  the  moment  coefficient 

caused  by  the  control  surface  deflection  about  either  the  l,m  ox  n  axis.  M  is  the  Mach 
number,  a  is  the  angle  of  attack  and  ^  is  the  control  surface  deflection  vector  caused  by 
the  eight  individual  vehicle  control  surfaces. 

The  first  task  of  the  control  allocator  determines  the  optimum  effector 
displacement  that  results  in  minimum  control  deficiency.  The  full  description  is 
contained  in  reference  22  and  repeated  here  for  completeness.  The  full  set  of  control 
surface  deflections  for  the  X-33  can  be  represented  by  the  eight-element  effector 
displacement  vector  previously  identified  in  equation  (1.1)  .  This  vector  is  subject  to  the 
limits  of  equations  (1.2)  and  (1.3)  or  a  subset  thereof  for  control  surface  failures.  The 
control  allocation  problem  uses  an  iterative  algorithm  that  queries  the  full  6-DOF  X-33 
aerodynamic  database.  The  piecewise  linear  constrained  control  allocation  problem  is 


MirigJ^  = 


[M,a,d)  +  {M,a,d)  +  C„^,  {M,a,d)  +  (M,«) 


subject  to 


Notice  that  the  upper  and  lower  bounds  on  d  are  dependant  upon  the  operating 
mode  of  the  vehicle.  A  vehicle  operating  under  nominal  conditions  has  the  full  range  of 
deflections  available  while  a  vehicle  operating  with  a  control  surface  failure  has  a  subset 
of  the  nominal  displacement  vector’s  range. 

A  nonzero  value  of  means  there  is  insufficient  control  power  to  rotationally 

trim  the  vehicle  (it  is  trim  deficient),  and  the  vehicle  will  depart  controlled  flight.  The 
trim  deficient  regions  in  Mach-alpha  space  can  be  mapped  and  identified  as  “keep-out” 

zones  for  the  outer  loop  trajectory  algorithm.  If  the  control  allocation  algorithm  is  able  to 
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rotationally  trim  the  vehicle,  any  excess  control  power  can  be  used  to  maneuver  the 
vehicle.  The  control  allocator  employs  a  second  control  allocation  optimization 
algorithm  to  minimize  the  deviation  from  a  minimum  trim-drag  condition.  This 
condition  is  represented  by  an  eight-element  preference  vector  set  to  zero.  The  second 
optimization  is 


MingJ^ 


subject  to  ^  ^  ^ 

C,AM,a,S)\  =  <i 


In  the  case  where  the  vehicle  is  operating  with  excess  control  power,  the  effector 
displacement  vector  is  the  output  of  the  second  optimization  routine,  otherwise  ^ 
minimizes  the  trim  deficiency.  The  effector  displacement  vector  that  is  the  solution  of 
the  control  allocator  is  represented  as  ^  to  delineate  it  from  the  infinite  combination  of 
control  surface  deflections  that  are  available  to  achieve  the  same  desired  control  but  are 
not  optimized  for  either  trim  deficiency  or  minimum  trim-drag. 

Given  S*  [a,  M)  (from  the  control  allocation  algorithm),  the  X-33  aerodynamic 

database  is  queried  over  the  Mach-alpha  grid  to  generate  the  wing-body  and  trim  lift  and 
drag  coefficients.  These  coefficients  when  summed  together  represent  the  total  lift  and 
drag  coefficients.  The  use  of  the  control  allocator  solution  effectively  decouples  the 
trajectory  generation  algorithm  from  the  inner  loop  and  allows  the  lift  and  drag 
coefficients  to  be  functions  of  angle  of  attack  and  Mach  number.  Effector  displacement 
remains  an  implicit  argument  and  allows  for  the  full  6-DOF  model’s  behavior  to  be 
captured  by  a  lower  order  model.  The  outer  loop  equation  (1.12)  is  elaborated  using  the 
lift  and  drag  coefficients  as  the  sum  of  the  wing-body  coefficients  plus  the  effector- 
dependent  trim  coefficients  using  the  control  allocator  solution. 
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Q(M,a)  =  Q  (M,a)  +  Q^,(M,«,^’  {M,a)) 
Co{M,a)  =  Co{M,a)  +  Cj,^(M,a,S*  {M,a)) 


(1.20) 


The  relationships  in  (1.20)  hold  for  a  vehicle  in  symmetric  flight  and  can  be 
stored  in  an  aerodynamic  lookup  table  for  each  Mach  and  angle  of  attack  combination 
subject  to  the  limits  on  the  effector  displacement  vector  dictated  by  the  presence  of 
nominal  or  failure  conditions.  Each  failure  condition  requires  a  separate  run  of  the 
control  allocator.  A  block  model  implementation  of  the  control  allocator  is  depicted  in 
Figure  8 


Aerodata  for 


^optimal 


Figure  8  AFRF  Control  Allocator 

3,  Trajectory  Replanning  Algorithm 

The  vehicle  is  assumed  to  operate  in  the  nominal  condition  until  a  failure  has  been 
identified  through  some  online  fault  detection  or  identification  algorithm  at  which  point 
the  lift  and  drag  table  corresponding  to  the  failed  vehicle  configuration  is  generated.  An 
approximate  block  diagram  of  the  guidance  loop  is  shown  in  Figure  9 
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Figure  9  Trajectory  Replanning 

4.  Failure  Mode 

Once  a  failure  has  been  identified,  the  PLP  control  allocator  generates  a  new  set 
of  control  deflections  for  the  failed  vehicle  and  uses  the  aerodynamic  database  to 
determine  the  new  trim  lift  and  drag  contributions.  In  the  process  of  computing  the 
control  deflections  for  the  failed  vehicles,  points  in  the  Mach-alpha  envelope  for  which 
the  vehicle  cannot  be  trimmed  are  identified.  These  trim  deficiency  regions  must  be 
avoided  and  are  thus  used  as  path  constraints  for  the  trajectory  replanning.  The  trajectory 
aerodynamic  data  table  is  changed  to  one  appropriate  for  the  vehicle  in  the  current  failed 
condition  and  a  new  trajectory  is  generated  for  the  new  aerodynamics  and  path 
constraints. 

5,  Failure  Mode  Effects 

The  trim  deficiency  path  constraint  on  the  available  range  of  angle  of  attacks  for  a 
given  Mach  number  requires  that  the  vehicle  must  be  rotationally  trimmable  (assuming 
no  side  slip)  for  a  given  angle  of  attack.  When  the  vehicle  operates  in  a  nominal 
configuration  it  is  capable  of  trimming  over  a  full  range  of  angle  of  attack,  however, 
when  the  vehicle  experiences  control  degradation,  its  available  range  of  angles  of  attack 
may  be  severely  reduced.  Doman,  et  al.  [ref  19]  give  full  treatment  of  this  subject.  Figure 
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10  depicts  the  available  values  of  angle  of  attack  that  will  allow  for  a  rotationally 
trimmed  vehicle  experiencing  no  failures.  Notice  that  because  the  maximum  trim 
deficiency  value  is  on  the  order  of  machine  precision  so  that  the  entire  range  of 
combinations  of  angles  of  attack  and  Mach  numbers  is  considered  to  be  available  in  the 
nominal  (unfailed)  case. 


-10  0  10  20  30  40  50 
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Figure  10  X-33  Trimmable  Values  for  Angle  of  Attack.(from  Ref  19.  ) 

Figure  1 1  shows  an  example  of  the  constrained  Mach-alpha  envelope  available  to 
the  X-33  following  a  dual  body  flap  failure  (stuck  at  26  degrees  up).  Notice  that  the 
envelope  describing  allowable  angles  of  attack  (in  dark  blue)  previously  depicted  in 
Figure  10  is  now  severely  reduced.  Angles  of  attack  other  than  those  contained  within 
the  dark  blue  region  will  place  the  vehicle  in  an  attitude  where  it  is  not  trimmable  in  pitch 
(pitch  deficient)  and  therefore  uncontrollable.  This  establishes  a  restricted  region  that 
feasible  trajectories  must  avoid.  In  other  words,  the  trajectory  must  have  a  clear  path  in 
Mach-alpha  space  from  the  initial  conditions  to  the  final  conditions.  A  clear  path  is 
defined  as  a  path  from  one  point  defined  by  the  flight  condition  (Mach-alpha)  to  another 
point  (Mach-alpha)  with  no  intervening  regions  of  unacceptable  trim  deficiency.  Should 
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a  path  cross  into  or  through  a  region  of  pitch  deficiency,  the  trajectory  becomes  infeasible 
because  the  vehicle  is  no  longer  rotationally  trimmable  and  therefore  uncontrollable. 

Pitch  Deficiency  Contour 
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Figure  1 1  X-33  Trimmable  Values  of  Angle  of  Attack  following  a  Dual  Stuck  Flap 

Failure  (failed  at  26  degrees  up),  (from  Ref  19). 

A  feasible  solution,  if  it  exists,  must  be  found  using  these  new  path  constraints. 
The  allowable  range  of  angle-of-attack  can  be  determined  in  a  manner  similar  to  the 
previous  method  used  to  generate  the  lift  and  drag  tables  using  the  AFRL  X-33 
aerodynamic  database  and  control  allocator.  The  result  is  a  dataset  that  accurately 
describes  the  trim  deficiency  as  a  function  of  Mach  number  and  angle  of  attack  (the 
vehicle’s  flight  condition).  This  dataset  is  the  output  of  the  AFRL  trim  deficiency 
algorithm  which  also  determines  the  new  coefficients  of  lift  and  drag  for  the  failed 
condition  using  the  new  optimized  effector  displacement  vector. 

Some  failures  will  have  very  restrictive  path  limits  for  feasible  trajectories  while 
others  may  have  no  clear  path.  Figure  12  shows  a  failure  where  the  vehicle  is  free  to 
operate  in  a  high  Mach  flight  condition  but  cannot  safely  land.  Figure  13  depicts  a 
failure  where  the  vehicle  is  limited  to  slower  speeds  and  alpha  ranges  to  safely  land. 
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Note  that  the  previous  diseussion  does  not  consider  any  other  limitations  such  as  normal 
force,  heating  or  dynamic  pressure  on  the  trajectory.  These  constraints  will  further 
reduce  the  area  available  to  a  clear  path  for  a  feasible  trajectory.  In  some  cases,  the 
additional  constraints  will  completely  eliminate  a  clear  path  and  render  trajectory 
infeasible.  It  is  easier  to  visualize  the  optimal  trajectory  with  constraints  as  a  region 
through  Mach-alpha  space  for  which  the  vehicle  must  be  rotationally  trimmable.  If  a 
failure  impinges  on  this  corridor  and  prevents  a  clear  path  through  to  landing,  then  there 
are  no  feasible  trajectories  available. 
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Figure  12  X-33  Pitch  Deficiency  Following  a  Stuck  Right  Rudder  (15  degrees) 
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Figure  13  X-33  Pitch  Deficiency  Following  a  Stuck  Flap  (at  -15  degrees) 

6,  Modeling  Trim  Deficiency 

As  the  previous  graphs  of  failure  modes  show,  the  trim  deficiency  can  be  modeled 
as  a  function  of  two  variables; 

d  =  f{M,a)  (1.21) 

The  function,  f,  can  be  determined  from  various  fits  to  the  trim  deficiency  data. 
This  thesis  uses  the  Table  Curve  3D  software,  version  4.0  by  Systat  Software 
{www.systat.com)  to  generate  MATLAB  function  codes  that  approximate  the  data 
generated  by  the  trim  deficiency  routines.  Table  Curve  3D  is  very  robust  and  can  apply 
and  analyze  the  fits  of  over  9000  equation  types  to  the  fitted  data  in  short  order.  An 
example  of  the  raw  data  from  the  trim  deficiency  code  imported  into  Table  Curve  3D  is 
shown  in  Figure  14 
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Trim  Deficiency  in  Mach-Alpha  Space 
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Figure  14  Table  Curve  3D  Trim  Deficiency  Data  for  a  Dual  Stuck  Flap  Failure  (26 

degrees) 

A  10*  order  Chebyshev  polynomial  produces  a  very  good  data  fit.  The  function 
generated  by  the  Table  Curve  3D  code  was  used  in  the  path  function  file  of  the  DIDO 
code.  The  path  function  file  combined  with  appropriate  bounds  on  the  path  prevents  the 
selection  of  combinations  of  Mach  and  angle  of  attack  that  exceed  a  given  trim 
deficiency.  An  important  selection  criterion  for  the  fitted  curve  is  the  function’s  behavior 
at  the  boundaries  of  the  Mach-alpha  envelope.  Some  fitting  functions  that  showed  good 
residuals  for  the  data  were  rejected  out  of  hand  because  the  low  Mach  and  angle  of  attack 
data  fits  were  very  poor.  It  is  more  important  that  the  fitting  function  models  the  data 
within  the  feasible  ranges  of  trim  deficiency  than  for  ranges  where  feasible  trajectories  do 
not  exist.  Mismatched  data  for  extreme  trim  deficient  regions  can  be  tolerated  but  it  is 
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essential  that  the  fit  near  the  transition  from  zero  trim  deficieney  to  positive  trim 
defieieney  be  as  good  as  possible. 


Trim  Deficiency  in  Mach-Alpha  Space 
Rank  1  Eqn  449  Chebyshev  X,LnY  Bivariate  Polynomial  Order  10 
r2=0.99796074  DF  Adj  r2=0.99790679  FitStdErr=0.0008263599  Fstat=  18791 .943 
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Figure  15  Trim  Defieieney  Fitted  with  10th  order  Chebyshev  Polynomial 


The  trim  defieieney  ean  also  be  modeled  as  a  pure  table  lookup  whieh  queries  the 
tablular  data  from  the  trim  deficieney  algorithm.  This  approach  is  more  accurate  than  the 
fitted  function  but  was  omitted  from  this  study  because  of  the  inherent  execution 
penalties  resulting  from  a  3-dimensional  table  lookup.  The  feasibility  of  following  a  trim 
deficiency  path  constraint  is  shown  by  projecting  the  trajectory  into  Mach-alpha  space 
along  with  the  trim-deficiency  contour  plot  and  verifying  that  the  trajectory  does  not 
violate  the  path  constraints.  The  feasibility  can  also  be  shown  by  inspection  of  the 
numerical  trim  deficiency  data  for  the  trajectory. 
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III.  THE  OPTIMAL  CONTROL  PROBLEM 


A.  OPTIMAL  CONTROL  THEORY 

The  following  sections  summarize  the  general  optimal  control  problem 
formulation  and  the  two  known  methods  of  generating  solutions. 

1.  Preliminaries 

Nearly  any  trajectory  generation  problem  can  be  posed  as  an  optimal  control 
problem  and  the  problem  considered  by  this  thesis  is  no  exception.  Given  a  system  with 
dynamic  constraints 

i  =  /(x,M,r)  (1.22) 

where 

/  :=  vector  of  functions  that  describe  the  dynamics  of  the  system 

X  :=  vector  of  states  that  fully  describe  the  system  at  any  r 

u  :=  vectors  of  control  variables 

T  :=  an  independent  variable  (usually  time) 


The  system  is  subjected  to  additional  constraints  to  include  the  path  constraint: 
hj  <  h[x,u,t)<  h^,  (1-23) 

where 

h  :=  vector  of  functions  describing  the  path  constraints 

hj  :=  vector  of  lower  path  bounds 

:=  vector  of  upper  path  bounds 
and  boundary  conditions 

5  ^e(x(ro),x(r^),ro,r^)<^  (1.24) 


where 
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vector  of  initial  boundary  conditions  (at  Tq  ) 


e{x{T^),T^)  := 

e{^x[T  ,T  :=  vector  of  final  boundary  conditions  (at  ) 
e^  :=  vector  of  lower  bounds 
:=  vector  of  upper  bounds 

The  system  also  has  bounds  on  the  control  and  state  variables  represented  by 


(1.25) 

^  <M(r)<^ 

(1.26) 

where 


x{t)  := 

state  vector  at  any  instance  of  r 

M(r)  := 

control  vector  at  any  instance  of  r 

:= 

vectors  of  lower  and  upper  state  bounds 

:= 

vectors  of  lower  and  upper  control  bounds 

The  bounds  on  the  state  and  control  vectors  can  be  thought  of  connecting  the 
physical  system  being  modeled  to  the  mathematical  model  used  to  solve  the  optimal 
control  problem.  For  example,  if  the  control  is  acceleration,  no  physically  realizable 
device  produces  infinite  acceleration  and  we  therefore  may  choose  to  limit  the  control  to 
a  sensible  value.  Similarly,  if  one  element  of  the  state  vector  is  velocity,  we  may  choose 
to  limit  the  velocity  to  one  that  is  within  the  design  limits  of  the  physical  system. 
Imposing  constraints  on  the  state  and  control  vectors  may  not  be  necessary  in  all  cases 
but  the  resulting  optimal  solution  of  unbounded  states  or  controls  needs  to  be  given  a 
hard  reality  check  when  realizing  a  physical  implementation  of  the  optimal  solution. 

An  optimal  control  problem  seeks  to  determine  the  solution  of  the  preceding 
dynamical  system  with  its  bounds  and  constraints  to  minimize  a  given  performance 
index.  The  performance  index  is  mathematically  defined  below: 
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where 


y(x(r),M(r),ro,r^)  =  £'(x(rj,x(?^.),r„,r^.)  + j7^(x(r),M(r),r)(ir  (1.27) 

^0 

E  :=  sealar  eost  funetion  evaluated  at  the  boundaries 

F  :=  sealar  eost  funetion  evaluated  over  the  entire  interval 

When  J  consists  of  only  E,  the  cost  function  is  in  Mayer  form.  When  J  consists  of  only 
F,  the  cost  function  is  in  Lagrange  form.  When  both  E  and  F  are  present,  the  cost 
function  is  in  Bolza  form. 

The  general  optimal  control  problem  is  fully  posed  in  the  following  manner; 

min„  7(x(r),M(r),r„,r^)  =  £'(x(rj,x(t^),r„,r^)  + jF(x(r),M(r),r)(ir 
subject  to  x  =  f  (x,  u,  r) 
hi  -  h[x,u,T)  < 

5  ^e(x(ro),x(r^),ro,r^)<^ 

Ui<u{t)<u^ 

The  formulation  of  the  system  as  represented  in  the  state  vector,  the  appropriate  choice  of 
the  dynamical  equations,  and  the  selections  of  the  constraints  and  bounds  forms  the  basis 
of  the  fully  posed  optimal  control  problem 

2.  Solution  Methodology 

Once  the  problem  has  been  properly  formulated  there  are  two  general  methods 
available  to  solve  the  optimal  control  problem,  direct  and  indirect  [ref.  27].  Indirect 
methods  tend  to  generate  fast  solution  times  and  greater  accuracy  but  are  considerable 
more  difficult  to  formulate  and  are  very  sensitive  to  the  initial  guess  [ref.  27]. 
Formulation  of  the  indirect  method  can  be  exceedingly  difficult  in  situations  where  the 
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dynamics  functions  are  not  pure  funetions  but  instead  rely  on  numerieal  data  such  as  the 
aerodynamie  lookup  tables  used  in  this  thesis. 

Direct  methods  reduce  the  optimal  eontrol  problem  to  a  single  large  Nonlinear 
Programming  (NLP)  problem.  The  strengths  of  the  indirect  methods  are  that  the 
formulation  is  significantly  easier  and  the  methods  are  relatively  insensitive  to  the  initial 
guess.  A  partieular  drawback  of  direct  methods  is  large  exeeution  times.  When  direet 
methods  are  used,  there  are  tradeoffs  between  exeeution  time  and  aeeuraey.  An  inerease 
in  the  size  of  the  NLP  program  does  not  necessarily  yield  a  proportional  inerease  in 
solution  aeeuraey  or  solution  time.  This  drawbaek  inhibits  real  time  optimal  control  using 
direct  methods  for  certain  problems  such  as  are  considered  in  this  thesis;  however,  as 
eomputing  power  continues  to  inerease,  the  potential  for  real-time  optimal  eontrol  will 
soon  be  realized. 

3.  DIDO 

This  thesis  uses  the  DIDO  numerieal  dynamie  optimization  software  developed 
by  I.M.  Ross  and  F.  Fahroo  of  the  Naval  Postgraduate  School.  DIDO  employs  a  direet 
Legendre  psuedospeetral  teehnique  that  uses  the  NLP  solver  SNOPT.  The  reader  is 
referred  to  References  12  through  17  for  a  full  description  of  the  Legendre  psuedospeetral 
method.  This  thesis  shows  another  application  of  DIDO  in  solving  optimal  trajectory 
problems. 

4,  Optimality 

Sinee  the  formulation  of  a  problem  using  a  direct  Legendre  psuedospeetral 
method  does  not  determine  the  adjoint  equations,  another  tool  must  be  employed  to  link 
the  solutions  generated  by  indireet  methods  and  direet  methods.  The  Coveetor  Mapping 
Theorem  (CMT)  provides  this  link  and  is  fully  explained  in  Reference  8.  The 
implieations  of  the  CMT  allow  a  determination  of  the  optimality  of  a  given  solution  by 
comparing  the  costates  of  the  solution  to  the  Kamsh-Kuhn-T ucker  (KKT)  multipliers. 
Every  time  a  state  or  control  becomes  bounded,  the  eorresponding  costate  should  follow 
the  KKT  conditions  as  if  they  were  the  KKT  multipliers.  Additional  evidence  of 
optimality  can  be  obtained  by  investigating  the  solution's  Hamiltonian  and  the  final  value 
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of  the  cost  function.  These  necessary  conditions  may  be  easily  obtained  from  the  DIDO 
solution  and  provide  insight  into  the  optimality  of  the  DIDO  generated  solution. 

5,  Feasibility 

Because  the  DIDO  solution  is  based  upon  the  controls  for  a  discretized  set  of 
points  determined  by  the  Legendre  polynomials  used  in  the  solution,  the  DIDO  solution 
may  not  possess  the  desired  accuracy  at  the  final  endpoint.  Since  the  accuracy  of  a  DIDO 
solution  is  a  function  of  the  number  of  nodes  evaluated  in  the  solution  space,  all  DIDO 
solutions  should  be  evaluated  for  feasibility.  This  is  usually  accomplished  by  comparing 
the  DIDO  solution  to  a  separate  solution  using  the  DIDO  controls  propagated  by  a  Runge 
Kutta  ODE  solver  such  as  MATLAB's  ODE45.  The  error  between  the  propagated  states 
and  the  DIDO-derived  states  gives  some  measure  of  confidence  of  the  accuracy  of  the 
DIDO  solution.  The  required  accuracy  is  determined  by  the  application,  for  example 
while  the  downrange  accuracy  of  a  20-node  DIDO  solution  reentry  trajectory  is  0.0001% 
of  the  propagated  downrange  distance,  it  still  may  be  inaccurate  enough  to  prevent  a 
feasible  landing  at  the  desired  field.  Eeedback  may  be  able  to  overcome  this  inaccuracy 
provided  the  solution  update  is  achieved  at  a  reasonable  rate.  Another  guard  against 
inaccuracies  in  a  trajectory  is  to  provide  additional  margins  when  defining  the  endpoint 
conditions. 

6.  Scaling  and  Balancing 

While  direct  methods  are  relatively  insensitive  to  the  initial  guess,  a  well- 
formulated  problem  can  be  completely  ruined  by  improper  scaling  and  balancing  of  the 
variables  and  constants  used  in  the  computation.  Ross  [ref  21]  notes  the  effect  of 
improper  scaling  on  the  DIDO  solution.  It  should  be  noted  that  inks  units  are  a  much 
better  choice  of  units  and  are  preferred  to  English  units  from  a  computational  viewpoint 
because  mks  units  are  direct  combinations  of  distance,  mass  and  time  base  units  and  thus 
easily  scaled  while  some  English  units  have  no  direct  equivalent. 

The  intent  of  scaling  and  balancing  is  to  avoid  the  problems  associated  with 
variables  whose  range  exceeds  that  of  numerical  precision  (about  16  significant  figures 
on  a  Pentium-class  PC).  A  proper  scaled  problem  scales  the  variables  to  a  range  of  0  to 
1 .  While  it  is  not  necessary  to  be  exactly  within  this  range,  it  is  clear  that  variables  of  the 
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same  order  of  magnitude  avoid  problems  associated  with  computing  ill-conditioned 
matrices  resulting  from  numerical  imprecision.  The  general  approach  to  scaling  the 
problem  variables  is  as  follows: 

The  fundamental  units  in  any  system  of  measure  traditionally  cover  mass,  time 
and  distance.  Other  units  such  as  area,  volume  and  velocity  are  derived  from  these 
fundamental  units.  Consider  the  initial  condition  for  the  third  element  of  the  state  vector, 
altitude.  The  initial  altitude  is  125,000  feet  while  the  final  altitude  is  500  feet. 
Additionally,  the  initial  downrange  distance  is  at  the  origin  of  the  coordinate  system 
while  the  expected  maximum  downrange  is  nearly  4,000,000  feet.  A  scaled  distance  unit 
is  created  by  applying  a  scale  factor  to  the  state  variables. 


X  = 

y  = 


X 

y 

DU 


z  = 


z 

DU 


Where  DU  represents  the  scaling  of  the  distance  unit.  To  scale  a  variable  that 
represents  velocity  notice  that 


,  .  distance 

velocity  = - 

time 

To  scale  velocity,  simply  divide  by  the  velocity  scaling  factor. 


V  = 


V 

VelUnit 


V 

'DU^ 

[Yuj 


=  v 


'  TU'" 
^~DUj 


Similarly,  the  scaling  factors  for  acceleration  and  force  are 


AcelUnits  = 


DU 
{TUf  ■ 


ForceUnits  =  MU 


DU 

Stu)" 


J 
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Arbitrarily  choosing  DU=100,000  ft  and  TU=  50  seconds  for  this  problem  seales 
the  state  variables  as  shown: 


X. 

K. 

Xj- 

Vf 

Zf 


0 

[a- 

“0“ 

0 

A 

0 

200,000/t 

f, 

1 

8000 ft!  see 

\Ji\ 

2 

4,000,000/r 

A 

~20 

0 

A 

0 

500/t 

A 

0.0025 

250/t/see 

vA 

0.0625 

After  sealing,  the  range  of  distanees  falls  between  0  and  twenty  for  down  range 
and  between  0  and  2  for  veloeity.  The  problem  should  also  be  balaneed,  that  is  all  path 
eonstraints  should  be  as  elose  to  the  same  order  of  magnitude.  While  this  may  not  be 
possible,  every  attempt  should  be  made  to  balanee  the  problem  as  well  as  to  seale  it. 

B,  REENTRY  VIEWED  AS  AN  OPTIMAL  CONTROL  PROBLEM 

After  sealing  and  balaneing,  the  atmospherie  entry  (or  reentry)  ean  be  formulated 
as  an  optimal  eontrol  problem. 

The  reentry  problem  eonsiders  a  vehiele  entering  the  atmosphere  from  an  initial 
position,  presumably  near  LEO  and  seeks  to  generate  an  optimal  trajeetory  that  brings  it 
to  a  final  state  speeified  by  an  endpoint  eondition.  The  trajeetory  may  not  violate  the 
path  eonstraints  whieh  limit  heating  rate,  normal  foree,  and  axial  foree  (dynamie) 
pressure.  The  optimal  trajeetory  seeks  to  minimize  a  performanee  index  whieh  is 
generally  a  funetion  of  the  state  of  the  vehiele.  Typieal  performanee  indiees  for  reentry 
are  maximizing  downrange,  maximizing  eross  range,  minimizing  eontrol  effort,  or  some 
weighted  eombination  of  these  indiees. 
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The  trajectory  is  also  subject  to  dynamic  constraints  which  describe  the  dynamics 
of  the  system.  As  discussed  in  the  previous  section,  this  thesis  uses  a  simple  dynamic 
model  that  retains  the  benefits  of  a  higher  fidelity  model  by  including  the  effects  of  trim 
lift  and  drag  and  using  table  lookups  into  a  high-fidelity  aerodynamic  database. 

For  a  typical  trajectory  that  seeks  to  determine  the  maximum  downrange  value  for 
an  entering  vehicle,  the  cost  function  is  simply  a  Mayer  cost: 

j(x(t),M(t),tofi/)  =  -X/ 

Time  is  the  independent  variable  for  this  problem.  A  Lagrange  component  of  the 
cost  can  be  included  as  desired.  The  initial  boundary  conditions  for  a  representative 
mission  are 


X.] 

“0 

L- 

0 

125,000/t 

Mach  8 

A 

-1.3  deg 

A 

0 

a, 

0 

0 

The  endpoint  equation  does  not  need  to  fully  describe  the  state  vector,  only  the  conditions 
that  are  necessary  for  a  safe  transition  and  landing.  The  manifold  can  be  simple  such  as: 


The  vector  above  does  not  specify  the  final  angles  of  attack  or  bank  and  it  is  possible  that 
the  vehicle  could  arrive  at  the  ending  manifold  in  an  unrecoverable  position  (high  angles 
of  attack  and  bank  along  with  large  vertical  speed).  To  prevent  this  from  occurring,  it  is 
good  practice  to  fully  specify  these  variables  as  well  as  the  associated  acceptable  range  of 
values  at  the  end  point.  Initial  trajectories  computed  appear  to  be  feasible  trajectories  but 


)  _ 

“500  ft 

_^f\ 
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upon  closer  inspection  met  the  end  point  manifold  with  vertical  speeds  exceeding  20,000 
feet  per  minute.  The  endpoint  manifold  was  then  modified  to  an  endpoint  set  limited  to 
vertical  speeds  between  -1500  ft/min  <  z  <  500  ft/min  . 

The  path  constraints  are  derived  primarily  from  the  design  and  mission  considerations  of 
the  vehicle. 
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normal  force(t) 
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1 
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There  are  limits  placed  on  the  state  vector  and  controls 
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40 
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The  limits  that  have  values  of  plus  and  minus  infinity  represent  states  without 
hard  limits  and  are  not  truly  constraints  upon  the  state  vector.  The  constraint  on  the 
altitude  component  of  the  state  vector  is  imposed  to  prevent  premature  contact  with  the 
ground. 
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IV.  APPLICATIONS 


A.  FOOTPRINT  DETERMINATION 

The  first  step  in  determining  the  vehiele  performanee  begins  with  investigating 
the  nominal  (unfailed)  ease.  A  rough  vehiele  footprint  ean  be  established  by  ehanging 
the  eost  funetions  to  inelude  various  eombinations  of  the  form 

M  =  Ax(^end)  + Byi^end) 

Where  A  +  B  =  1  and  0<y4<l;0<5<l. 

x[end)  and  y[end)  are  the  final  values  of  downrange  and  erossrange. 

The  eonstants  for  eost  funetion  used  to  determine  the  maximum  downrange  are 
simply  =7  it  follows,  therefore  that  B=0.  Note  that  when  y4=-7,  minimum  downrange  is 
implied.  The  eonstants  for  eost  funetion  used  to  determine  the  maximum  erossranges  are 
simply  B=1  it  follows,  therefore  that  A=0.  Note  that  when  B=-l,  minimum  erossrange  is 
implied,  whieh  for  the  nominal  ease  should  be  symmetrie  about  the  x  axis. 

B,  LANDING  AT  A  SPECIFIC  FIELD  WITHIN  THE  FOOTPRINT 

Onee  the  rough  limits  of  the  footprint  have  been  established,  the  vehiele  should  be 
able  to  land  at  any  field  within  the  footprint.  In  order  to  reeover  the  vehiele  at  a  speeifie 
field,  the  eonditions  at  the  landing  threshold  that  lead  to  a  “sueeessful”  landing  must  be 
strietly  defined.  Original  runs  that  speeified  a  small  subset  of  the  state  veetor  at  the 
endpoint  led  to  unsatisfaetory  eonditions  at  the  landing  threshold.  Consider  the  endpoint 
eonditions  speeified  by 

e(xf,t^)  = 

Optimal  trajeetories  ean  be  generated  whieh  meet  these  endpoint  eonditions; 
however,  the  vehiele  state  may  be  sueh  that  a  sueeessful  landing  is  not  possible  given  the 
praetieal  limitations  of  the  vehiele.  Original  trajeetories  whieh  used  the  previous 
endpoint  eondition  appeared  reasonable  until  the  vertieal  speed  at  the  landing  threshold 
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was  investigated.  Some  trajeetories  hit  the  landing  threshold  with  vertieal  speeds  in 
exeess  of  4000  feet/see.  Clearly,  this  is  an  unreeoverable  attitude-  the  vehiele  passes 
through  the  endpoint  on  its  way  to  its  final  resting  plaee  at  the  bottom  of  a  very  deep 
erater  instead  of  the  end  of  the  runway.  Another  eonsideration  is  angle  of  bank.  Hitting 
the  terminal  eonditions  with  strong  left  drift  and  extreme  right  bank  angle  ean  be 
disastrous.  Limits  on  the  bank  angle  are  elearly  needed.  The  endpoint  eonditions  that 
speeify  a  landable  attitude  have  now  grown  to 
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Notiee  that  angle  of  attaek  is  not  speeified  in  the  endpoint  eonditions  and  there  is 
a  ehanee  that  the  vehiele  eould  be  near  stall.  The  endpoint  eonditions  as  stated  above  are 
insuffieient  to  hit  a  speeified  point  along  the  trajeetory.  The  landing  field  is  speeified  by 
Cartesian  eoordinates  x,y,z  as  well  as  a  runway  heading,  speeified  by  the  azimuth  angle, 
P .  Notiee  that  in  the  flat  earth  model,  the  altitude  of  the  field  is  assumed  to  be  zero  but 
is  not  neeessarily  zero.  Endpoint  eonditions  for  landing  at  fields  with  nonzero  elevations 
are  simply  aeeommodated  as 

~  ^ field  ^ojfset 

The  new  endpoint  eonditions  for  landing  at  a  speeified  field  in  a  landable  vehiele 
attitude  are  now  speeified  as 

=  T/  Vf  Pf 


The  runway  heading  for  the  seleeted  landing  field  ean  have  signifieant  effeet  on 
the  trajeetory  beeause  reversing  direetion  requires  eonsiderable  energy  loss  during  the 
trajeetory.  It  is  entirely  possible  that  a  footprint  generated  with  loosely  speeified 
endpoint  eonditions  does  not  eontain  interior  points  for  feasible  trajeetories  when  the 
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endpoint  conditions  are  more  fully  specified  (such  as  a  runway  heading  180  degree  from 
the  original  vehicle  azimuth). 
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V.  RESULTS 


A.  DISCUSSION 

The  following  sections  discuss  the  footprint  and  individual  trajectory  analysis  for 
several  cost  functions.  Unless  otherwise  specified,  the  trajectories  for  both  failed  and 
nominal  vehicles  use  the  same  initial  conditions  of 

e(x,,t,)  =  [x,  y,  z.  V,  y,  a, 

With  the  following  values 

e(x.,t.)  =  [0  ft  Oft  125  kft  Mach  8  -1.3°  0°  0°  0°f 

These  initial  conditions  correspond  to  a  point  on  the  entry  trajectory  after  the  peak 
heating  and  dynamic  pressure  and  are  suitable  for  failure  as  the  vehicle  can  begin  to  use 
its  control  surfaces  to  control  the  trajectory.  The  initial  condition  for  the  angle  of  attack 
may  be  too  shallow  based  on  previous  work,  but  it  is  of  little  consequence  for  studying 
the  effects  of  failure  on  trajectory  replanning  and  retargeting. 


The  endpoint  conditions  for  both  trajectories  are  of  the  form 

=  Zf  Zf 

With  the  following  values 

=  [free  free  free  -500 ±1000  ft/ min  Mach  0.15 ± 0.003  free  0±20°]^ 

Both  vehicles  observe  a  normal  force  constraint  of  3±/-  6  g’s.  The  vehicle  experiencing 
a  control  surface  failure  must  observe  the  trim  deficiency  path  constraint  while  the 
nominal  vehicle  is  unconstrained  by  trim  deficiency  because  as  shown  in  Figure  10  ,  the 
vehicle  has  adequate  control  power  to  trim  over  the  entire  flight  envelope. 

Once  the  rough  footprint  boundaries  are  established  for  the  vehicle  experiencing 
control  surface  failures,  the  algorithm’s  ability  to  generate  trajectories  to  a  specific 
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landing  field  and  its  associated  runway  is  demonstrated.  For  this  case,  the  Cartesian 
coordinates  and  azimuth  are  no  longer  free  and  are  replaced  with  the  field  xy  coordinates, 
an  altitude  corresponding  to  500  feet  above  the  ground  and  a  corresponding  runway 
azimuth. 

B,  FOOTPRINT  ANALYSIS  FOR  NOMINAL  AND  STUCK  FLAP 

TRAJECTORIES 

An  approximate  footprint  for  the  nominal  and  stuck- flap  trajectories  can  be 
determined  by  choosing  the  cost  functions  corresponding  to  minimum  and  maximum 
down  range  distance  as  well  as  minimum  and  maximum  cross  range  distance.  The  four 
points  form  an  approximate  boundary  of  the  footprint  and  feasible  trajectories  to  the 
interior  points  of  the  polygon  formed  by  these  boundary  points  should  exist.  Note  that 
feasible  trajectories  to  exterior  points  are  not  ruled  out  but  should  be  proven  for  each 
case.  Additional  boundary  points  can  be  determined  by  varying  the  A  and  B  terms  in  the 
previous  equation  M  =  Ax  [  end )  + By  [  end ) . 

1,  Aerodynamic  Considerations 

The  introduction  of  a  control  surface  failure  can  drastically  change  the 
aerodynamics  of  the  vehicle.  The  coefficients  of  lift  and  drag  for  the  vehicle 
experiencing  a  stuck  flap  have  changed  significantly.  The  constraints  on  angle  of  attack 
resulting  from  the  trim  deficiency  requirements  produce  a  vehicle  which  may  look  the 
same  but  flies  considerably  different.  The  footprints  for  the  failed  and  normal  trajectories 
should  appear  different  because  of  the  different  aerodynamic  configurations  and  Mach- 
alpha  combinations. 

Despite  the  aerodynamic  changes  of  the  failed  vehicle,  its  performance  should  be 
a  subset  of  the  nominal  vehicle’s  performance.  This  is  a  valid  assumption  because  the 
control  surface  failures  encountered  in  flight  limit  the  control  surface  deflections  to  a 
subset  of  their  nominal  values.  In  other  words,  the  nominal  vehicle  can  fly  the  same 
profile  as  a  vehicle  with  control  surface  failures  simply  by  voluntary  limiting  the  unfailed 
control  surface  deflections  to  those  of  the  vehicle  with  the  failures.  The  only  difference 
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becomes  semantic;  the  nominal  vehicle  “simulates”  the  failure  while  the  vehicle  with  the 
actual  failure  has  no  choice  but  to  fly  (if  it  can)  with  the  current  failed  conditions 

2,  Footprint  Observations  and  Conclusions 

Figure  16  shows  the  footprints  for  four  selected  trajectories:  nominal,  stuck  flap 
failure  with  loosely  specified  endpoint  conditions,  stuck  flap  failure  with  tightly  specified 
endpoint  conditions  and  stuck  flap  failure  with  tightly  specified  endpoint  conditions  and 
trim  deficiency  path  constraint.  Several  things  are  readily  apparent: 

a.  The  footprints  for  vehicles  with  failures  have  better  minimum  downrange 
performance  than  that  of  the  nominal  vehicle.  This  was  surprising  but  can  be  explained 
by  the  control  allocation  algorithm.  Recall  that  for  a  given  Mach-alpha  pair,  the  control 
allocator  computes  a  vector  of  control  deflections  that  simultaneously  balance  the 
moments  produced  by  the  wing-body  and  minimize  the  1-norm  of  the  control  deflection 
vector.  Changing  the  preference  vector  to  maximize  trim-drag  or  control  deflections 
within  the  control  allocator  can  solve  this  particular  issue;  however,  the  differences  in  the 
footprints  are  indicative  of  greater  problems  associated  with  having  the  outer  loop 
decoupled  from  the  inner  loop. 

b.  There  is  great  variance  between  the  three  footprints  of  the  failed  vehicles. 
The  largest  footprint  for  the  failed  vehicles  represents  the  traditional  footprint  for  a  failed 
vehicle  as  found  in  the  literature.  One  such  approach  to  determining  this  footprint  uses 
Vinh’s  [Ref  23]  equations.  The  footprint  uses  loosely  specified  end  point  conditions 
(altitude  and  airspeed),  and  does  not  prevent  the  angle  of  attack-Mach  combination  from 
entering  regions  of  the  flight  envelope  where  the  vehicle  cannot  achieve  moment 
equilibrium  due  to  a  lack  of  control  power.  As  such,  the  larger  footprint  depicts  where 
the  “pieces  fall”  since  in  the  early  moments  following  a  failure,  the  trajectory  would 
violate  the  trim  deficiency  constraint  and  the  vehicle  would  become  uncontrollable  and 
break  up. 

c.  The  effect  of  more  restrictive  endpoint  conditions  further  reduces  the 
footprint.  The  third  largest  footprint  represents  the  footprint  for  the  trajectory  for  the 
failed  vehicle  that  observes  the  trim  deficiency  path  constraint.  While  this  trajectory 

appears  to  be  sufficient,  further  inspection  reveals  that  the  endpoint  conditions  are  met 
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with  a  vertical  speed  of  nearly  20,000  feet  per  minute.  The  vehicle  meets  the  endpoint 
conditions  on  its  way  to  a  deep  crater  as  there  is  insufficient  remaining  altitude  to  slow 
the  rate  of  descent  to  an  acceptable  level.  The  smallest  footprint  represents  a  trajectory 
for  the  failed  vehicle  with  tightly  specified  endpoint  conditions  that  limit  the  vertical 
speed  and  angle  of  bank  to  levels  that  permit  landing.  This  footprint  for  a  “landable 
vehicle”  in  the  failed  condition  is  1%  of  the  area  of  the  largest  failed  footprint  and  shows 
that  trajectories  generated  using  energy  analysis  with  loosely  specified  endpoint 
conditions  and  without  trim  deficiency  path  constraints  result  in  overly  optimistic 
footprint  estimates. 


Footprint  Comparisons  of  Various  Trajectories 


Figure  16  Comparison  of  Nominal  and  Stuck-Flap  (+26°)  Footprints 
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Footprint  Comparisons  of  Various  T rajectories 


Figure  17  Contraction  of  Failure  Stuck-Flap  (+26°)  Trajectory  Footprint 


The  difference  between  a  “landable”  and  “unlandable”  footprint  is  shown  in 
Figure  17  Notice  that  the  introduction  of  vertical  speed  into  the  endpoint  set  reduces  the 
footprint  noticeably.  Landable  footprints  result  from  trajectories  with  endpoint 
conditions  that  allow  a  safe  transition  from  the  landing  threshold  to  weight  on  wheels  and 
rollout  that  are  within  the  vehicle  and  payload  performance  parameters. 

3,  Maximum  Downrange 

The  following  trajectories  for  the  nominal  and  failed  vehicle  were  conducted  as 
part  of  the  rough  footprint  analysis.  The  performance  criterion  for  these  trajectories  was 
to  maximize  the  downrange  position  at  the  end  of  the  trajectory.  The  following  trajectory 
analysis  is  a  detailed  comparison  of  the  nominal  and  failed  vehicle  trajectories  for  each 
associated  cost  function. 
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a.  Mach-Altitude  Profile 

The  Mach-altitude  profiles  for  both  trajectories  are  shown  in  Figure  18  . 
Notice  that  the  nominal  trajectory  trades  excess  speed  for  altitude.  At  the  top  of  the  zoom 
climb,  the  vehicle  assumes  the  best  lift/drag  configuration  and  continues  to  the  surface. 
The  trajectory  for  the  vehicle  with  the  stuck  flap  is  prevented  from  using  the  same 
combination  of  angle  of  attacks  as  the  normal  because  of  the  imposed  trim  deficiency 
requirement.  The  trajectory  does  use  the  largest  allowable  angle  of  attack  but  this  is 
insufficient  to  produce  the  zoom  climb  observed  in  the  nominal  trajectory. 


Mach  number 


Figure  18  Maximum  Downrange  Mach- Altitude  Profile 
b.  Trajectory  in  Cartesian  Coordinates 

Figures  19-22  depict  the  Cartesian  coordinates  of  the  trajectory.  The  x,  y, 
and  z  coordinates  are  direct  results  of  the  controls  and  the  time  of  their  applications.  The 
only  performance  measure  for  these  trajectories  is  to  maximize  the  downrange  (x).  There 
is  no  limit  placed  on  crossrange  (y)  and  the  altitude  must  meet  the  endpoint  condition  of 
500  feet. 
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The  downranges  for  the  failed  and  nominal  trajectories  differ  by  nearly 
200  miles.  This  is  not  unreasonable  considering  that  the  symmetrical  flap  failure 
drastically  increases  the  drag  on  the  vehicle.  This  drag  when  combined  with  the  imposed 
trim  deficiency  path  constraint  considerably  slows  the  vehicle  and  prevents  it  from 
trading  altitude  for  velocity.  The  total  trajectory  time  of  flight  is  much  shorter. 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  19  Downrange  Profile  for  Maximum  Downrange  Trajectories 

The  crossrange  profiles  of  both  trajectories  exhibit  some  unexpected  left 
drift.  This  drift  is  insignificant  when  compared  with  the  total  downrange  distance  of  the 
trajectory.  The  cause  of  the  drift  is  unknown  and  may  be  numerical.  From  a  physical 
point  of  view,  the  aerodynamics  of  both  the  failed  and  nominal  vehicles  is  symmetrical  as 
is  the  model  used  for  the  trajectory.  Closer  examination  of  the  control  allocation  and 
aerodynamic  data  may  yield  answers.  One  obvious  difference  between  the  two 
trajectories  is  that  the  nominal  trajectory  drifts  slowly  to  the  left  while  the  vehicle  with 


49 


stuck  flaps  experiences  both  left  and  right  drift  along  the  trajeetory.  The  angle  of  bank 
should  be  inspeeted  to  eorrelate  the  drift  in  erossrange. 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  20  Cross  Range  Profiles  for  Maximum  Downrange  Trajectories 

The  altitude  profdes  have  been  diseussed  in  some  detail  previously,  but  an 
interesting  event  oeeurs  near  the  end  of  eaeh  trajeetory.  Both  trajeetories  overshoot  the 
endpoint  altitude  of  500  feet  and  appear  to  perform  identieal  maneuvers  to  meet  the  500 
foot  altitude  requirement  at  the  endpoint  condition.  This  event  was  surprising  not  only 
because  of  the  nearly  identical  maneuver,  but  also  beeause  the  trajectory  actually 
indicated  a  climb  in  the  last  moments  of  the  trajectory.  The  encroaehment  below  500 
foot  altitude  was  allowed  by  the  problem  formulation;  the  only  limit  plaeed  upon  the  z 
state  variable  was  that  it  had  to  be  nonnegative. 
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Altitude,  feet 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  21  Terminal  Maneuver  for  Max  Downrange  Trajectories 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  22  Altitude  Profiles  for  Max  Downrange  Trajectories 
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c.  Velocity  Profiles 

Figure  23  depicts  the  velocity  profiles  for  the  max  downrange  trajectories. 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  23  Velocity  Profiles  for  Max  Downrange  Trajectories 


d.  Flight  Path  and  Azimuth  Angle 

The  flight  path  angles  (FPA)  shown  in  Figure  24  for  these  two  trajectories 
differ  as  a  consequence  of  the  controls  available  to  the  individual  vehicles.  The  failed 
vehicle  exhibits  a  very  steep  descent  but  ends  the  trajectory  with  a  terminal  maneuver 
similar  to  the  nominal  trajectory.  The  similarity  of  the  terminal  maneuvers  is  a  direct 
result  of  the  endpoint  conditions  shared  by  the  two  trajectories. 
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Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


time,  seconds 

Figure  24  Flight  Path  Angle  Profiles  for  Max  Downrange  Trajectories 

There  is  no  limit  placed  upon  the  azimuth  angles  for  these  trajectories  and 
both  vehicles  wander  somewhat  in  azimuth.  The  average  azimuth  of  both  trajectories  in 
Figure  25  is  positive  and  accounts  for  the  left  bias  of  the  crossrange.  Notice  the  nominal 
trajectory  has  the  most  stable  azimuth  between  100  and  700  seconds.  The  stable  azimuth 
is  small  and  positive  and  accounts  for  the  steady  left  drift  found  in  the  nominal  trajectory. 
The  vehicle  experiencing  the  stuck  flap  is  far  less  stable  and  the  azimuth  varies  by  as 
much  as  1 1  degrees. 
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Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 
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Figure  25  Azimuth  Angle  Profile  for  Max  Downrange  Trajeetories 
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e. 


Vertical  Speed  Profile 


Nominal  Trajectory 


Trajectory  with  Flaps  Stuck  at  26  degrees 


Figure  26  Vertical  Speed  Profiles  for  Nominal  and  Failed  Trajectories  (Max 

Downrange) 

f.  Body  Frame  Normal  Force 

The  normal  force  experienced  by  the  vehicle  for  both  trajectories  are 
depicted  as  functions  of  time  and  distance.  As  expected,  the  vehicle  with  the  control 
surface  failure  has  a  more  varied  normal  force  profile  with  larger  force  excursions.  Both 
vehicles  have  increases  in  normal  force  at  the  terminal  maneuver  as  shown  in  Figure  27  . 
Additionally,  the  nominal  vehicle  is  “pulling  g’s”  during  the  zoom  climb  maneuver  early 
in  the  trajectory. 
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normal  force,  (in  gs)  c  normal  force,  (in  gs) 


Symmetrical  Stuck  Flap  Trajectory  80-nodes 


27  Normal  Force  ( - 1 )  vs.  Time  for  Vehiele  with  Stuck  Flaps 
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Figure  28 


Flaps 


Nominal  Trajectory  80-nodes 


Figure  29  Normal  Force  (  - 1 )  vs.  Time  for  Nominal  Vehicle 

Nominal  Trajectory  80-nodes 


Figure  30  Delta  Normal  Force(  n^-l)  vs.  Downrange  Distance  for  Nominal  Vehicle 
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g.  Controls 

While  the  problem  formulation  uses  pseudo  controls  to  induce  inertia  into 
the  dynamical  model,  the  physical  controls  of  the  vehicle  are  its  angles  of  attack  and 
bank.  The  angle  of  attack  profiles  clearly  demonstrate  the  effects  of  the  reconfigured 
aerodynamics  and  resulting  angle  of  attack  limits.  The  failed  vehicle  is  limited  to  a  much 
smaller  range  of  angles  of  attack  until  the  subsonic  flight  regime  at  which  point  it  can  use 
the  same  full  range  of  angles  of  attack  as  the  nominal  trajectory.  Both  trajectories  share 
the  same  terminal  maneuver.  The  terminal  maneuver  is  simply  a  climb  that  compensates 
for  undershooting  the  500  foot  final  altitude  condition.  This  maneuver  can  be  compared 
to  the  terminal  maneuver  shown  in  the  altitude  profile. 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  3  f  Angle  of  Attack  Profiles  for  Max  Downrange  Trajectories 

The  bank  angle  profiles  are  as  expected  and  compare  well  to  the 
crossrange  profiles  previously  discussed.  The  nominal  vehicle  is  relatively  steady  all  the 
way  to  the  endpoint  conditions  where  it  gets  a  little  chatter  while  the  vehicle  with  the 
stuck  flaps  has  bank  angle  excursions  as  large  as  27  degrees. 
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Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 
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Figure  32  Bank  Angle  Profiles  for  Max  Downrange  Trajectories 


h.  Pseudo  Controls 

The  pseudo  controls  depicted  are  not  the  true  controls,  but  rather  the 
control  rates.  The  profiles  show  that  while  the  majority  of  the  maneuvers  request  the 
controls  at  the  maximum  rate  available  (the  pseudo  controls  are  rate  limited),  not  every 
maneuver  requires  full  control  surface  rates.  The  points  where  the  pseudo  controls  are 
limited  provide  excellent  insight  into  the  optimality  of  the  solution  by  examination  of  the 
necessary  conditions. 
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Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  33  Angle  of  Attaek  Rate  Profiles  for  Maximum  Downrange  Trajeetories 


time,  seconds 


Figure  34  Rank  Angle  Rate  Profiles  for  Maximum  Downrange  Trajeetories 
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L  Trim  Deficiency 

The  restricted  region  in  Mach-alpha  space,  due  to  the  moment  trim 
deficiency  for  the  vehicle  with  stuck  flaps,  imposes  severe  trajectory  performance 
penalties  as  previously  shown.  Figure  35  is  meant  to  show  the  differences  between  the 
two  trajectories  in  Mach-alpha  space  as  well  as  to  provide  visual  verification  that  the 
trajectory  for  the  vehicle  experiencing  flap  failure  does  not  violate  the  imposed  trim 
deficiency  constraint  of  0.003  and  instead  closely  follows  the  contour  in  regions  where 
the  constraint  is  active. 


T rajectory  Trim  Deficiency 


Figure  35  Trim  Deficiency  Map  for  Max  Downrange  Trajectories 


j.  Optimality  Analysis 

The  optimality  of  the  two  trajectories  can  be  shown  by  investigating  the 
Hamiltonian  and  necessary  conditions.  The  Hamiltonian  for  these  trajectories  should  be 
zero.  Several  problems  inherent  with  the  trajectory  problem  as  posed  can  prevent  DIDO 
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from  identifying  a  locally  optimal  solution  and  instead  characterize  the  solution  as  “near 
optimal”  where  the  cost  may  be  within  99.999%  of  optimal.  The  two  major  problems 
which  lead  to  “near  optimal”  solutions  are 

1.  The  use  of  nonsmooth  data.  In  this  case,  the  aero  data  tables  are  not 
smooth  and  can  cause  some  difficulty  in  the  Jacobian  matrices  used  in  the  DIDO  code. 

2.  Insufficient  number  of  nodes  used  in  the  solution.  This  is  a  tradeoff 
between  execution  time  and  accuracy.  If  a  40-node  solution  yields  propagation  errors  of 
acceptable  magnitude  along  with  a  near-optimal  solution,  it  may  be  undesirable  to 
increase  the  execution  time  to  find  a  locally  optimal  solution,  if  it  exists. 

Both  nominal  and  stuck  flap  trajectories  used  an  80-node  solution  which 
resulted  in  a  near-optimal  solution  determination.  The  Hamiltonians  are  not  exactly  zero 
but  are  close.  The  stuck  flap  trajectory  Hamiltonian  is  not  as  smooth  as  expected  but 
investigation  of  the  necessary  conditions  shows  appropriate  switching  of  the  dual  for 
constrained  controls. 


80-node  Nominal  Trajectory,  Max  Down  Range 


Figure  36  Hamiltonian  for  Max  Downrange  Nominal  Trajectory 
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aoa  rate,  radians/sec  Hamiltonian 


Symmetrical  Stuck  Flap  Trajectory,  80-nodes 


Figure  37  Hamiltonian  for  Vehicle  with  (26  deg)  Stuck  Flaps 

Symmetrical  Stuck  Flap  Trajectory,  80-nodes 


Figure  38  Constrained  Pseudo  Controls  for  (26  deg)  Stuck  Flap  Trajectory 
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The  pseudo  controls  for  the  vehicles  experiencing  the  stuck  flaps  are 
generally  constrained  throughout  the  trajectory  as  shown  in  Figure  38  and  Figure  39  The 
rapid  switching  between  extreme  pseudo  control  rates  is  a  close  approximation  of  bang- 
bang  control.  The  necessary  conditions  for  the  pseudo  controls  are  satisfied  when  the 
duals  for  both  angle  of  attack  and  bank  angle  rates  are  positive  for  upper-bounded  control 
rates  and  negative  for  lower-bounded  control  rates  as  shown  in  Figure  40  and  Figure  41 


Symmetrical  Stuck  Flap  Trajectory,  80-nodes 


Figure  39  Upper  and  Lower  Bounded  Bank  Angle  Rate  for  (26  degree)  Stuck  Flap 

Trajectory 
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controls  dual 


Figure  40  Rank  Angle  Rate  and  Associated  Dual  for  (26  deg)  Stuck  Flap  Trajectory 


controls  dual 


controls 


Figure  41  Angle  of  Attack  Rate  and  Associated  Dual  for  (26  degree)  Stuck  Flap 

Trajectory 
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k.  Feasibility 

The  feasibility  of  the  trajectory  is  demonstrated  by  propagating  the  DIDO 
controls  using  an  R-K  algorithm  and  comparing  the  DIDO  and  propagated  state  vectors  at 
the  endpoint.  There  is  an  obvious  tradeoff -a  trajectory  that  seeks  to  place  the  vehicle  in 
a  landable  attitude  at  the  landing  threshold  of  a  runway  needs  sufficient  accuracy  at  the 
endpoint  conditions  that  does  not  place  the  vehicle  in  danger.  The  accuracy  of  an  80- 
node  solution  is  generally  sufficient.  The  errors  in  the  corresponding  components  of  the 
state  vector  are  shown  below. 

error  =  [436ft  -267ft  79  ft  11.5  deg  -10  deg  -21  deg  0.11  deg/sec  -0.63  deg/sec]’ 

A  graphical  depiction  of  the  DIDO  and  propagated  controls  is  shown  in 

Figure  42 


DIDO  and  Propagated  Controls  for  Stuck  Flap  T rajectory  -  80  nodes 


time,  secs 


Figure  42  DIDO  and  Propagated  Controls  (DIDO  controls  are  denoted  by  discrete 

points) 
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4,  Minimum  Downrange 

a.  Mach-Altitude  Profile 

The  Mach- altitude  profile  for  the  minimum  downrange  trajectory 
produced  some  unexpected  results.  The  nominal  and  stuck  flap  trajectories  closely 
approximate  each  other  until  approximately  60,000  feet  where  the  nominal  vehicle 
conducts  a  zoom  climb  that  approaches  stall  and  then  rapidly  descends  at  a  high  subsonic 
speed.  The  point  of  the  climb  corresponds  to  altitude  of  minimum  speed  of  sound  in 
Figure  7 


Mach  number 


Figure  43  Mach-altitude  Profiles  for  Minimum  Downrange  Trajectory 
b.  Trajectory  in  Cartesian  Coordinates 

The  profiles  in  Cartesian  coordinates  for  the  minimum  downrange 
trajectory  are  shown  in  Figure  44  .  The  stuck  flap  trajectory  shows  better  minimum 
downrange  performance  than  the  nominal  trajectory.  This  is  surprising  because  the 
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nominal  vehicle  has  the  full  available  range  of  control  effector  positions  available  while 
the  failed  vehicle  only  has  a  subset  of  the  available  control  effector  positions.  In  other 
words,  the  performance  of  the  nominal  vehicle  should  never  be  less  than  that  of  the 
vehicle  experiencing  a  control  surface  failure.  One  explanation  for  the  performance 
difference  lies  with  the  control  allocation  found  in  the  inner  loop  of  the  algorithm.  Recall 
that  the  control  effector  position  vector  for  the  requested  angle  of  attack  is  optimized  to 
minimize  pitching  moment  using  the  aero  dynamic  data  and  vehicle  flight  condition 
(determined  by  Mach  and  angle  of  attack)  as  arguments.  The  resulting  control  effector 
position  vectors  are  different  with  each  aerodynamic  configuration  for  the  same  requested 
angle  of  attack  with  accompanying  differences  in  the  lift  and  drag  forces  experienced  by 
the  vehicles.  The  consequence  is  nearly  10  miles  in  additional  downrange  distance  and 
100  additional  seconds  of  flight  time. 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  44  Downrange  Distance  Profiles  for  Minimum  Downrange  Trajectories 

The  crossrange  profile  of  the  vehicle  with  stuck  flaps  veers  nearly  a  mile 
off  centerline  before  attempting  to  return  to  centerline  while  the  nominal  trajectory 
remains  on  centerline  throughout  the  trajectory. 
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Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 
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Figure  45  Crossrange  Profiles  for  Minimum  Downrange  Trajeetories 


The  altitude  profile  for  the  nominal  trajeetory  shows  the  previously 
pointed  out  climb  at  50,000  feet.  Both  trajectories  exhibit  similar  terminal  maneuvers  as 
those  meeting  the  maximum  downrange  criterion.  The  vehicles  pass  through  500  feet 
and  initiate  climbs  to  meet  the  endpoint  conditions. 
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Altitude,  feet 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  46  Altitude  Profile  for  Minimum  Downrange  Trajeetories 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  47  Terminal  Altitude  Maneuver  for  Minimum  Downrange  Trajeetories 
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c.  Velocity  Profile 

The  vertical  speed  profiles  are  nearly  identical  until  the  transonic  region  at 
which  point  they  diverge.  This  divergence  is  a  consequence  of  the  nominal  trajectory’s 
zoom  climb  that  occurs  at  50,000  feet.  The  final  30  seconds’  velocity  profile  for  both 
trajectories  is  very  similar. 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  48  Velocity  Profiles  for  Minimum  Downrange  Trajectories 
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Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  49  Terminal  Velocity  Profile  for  Minimum  Downrange  Trajectories 


d.  Flight  Path  and  Azimuth  Angle 

The  terminal  flight  path  angle  profile  for  both  trajectories  is  very  similar. 
The  failed  vehicle’s  final  azimuth  performance  is  far  better  than  the  nominal  vehicle  with 
some  attempt  at  ground  track  reversal  in  evidence. 
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Azimuth  Angle,  degrees  Ciq’  Flight  Path  Angle,  degrees 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 
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Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 
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Figure  5 1  Azimuth  Angle  Profile  for  Min  Downrange  Trajectories 
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e.  Vertical  Speed  Profile 


Nominal  Trajectory 


Figure  52  Vertical  Speed  Profiles  (vs.  Downrange  Distance)  for  Minimum 

Downrange  Trajectories 
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Trajectory  with  Flaps  Stuck  at  26  degrees 


Figure  53  Vertical  Speed  Profiles  (vs.  time)  for  Nominal  and  Failed  Trajectories, 

Minimum  Downrange 
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normal  force,  (in  gs) 


f.  Body  Frame  Normal  Force 

Normal  Force  for  Stuck  Flap  Trajectory,  80  nodes 


Figure  54  Body  Frame  Normal  Foree  Stuek  Flap  Minimum  Downrange 

Trajeetory 
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Normal  Force  for  Nominal  Trajectory,  80  nodes 


Figure  55  Body  Frame  Normal  Foree  Nominal  Minimum  Downrange 

Trajectory 
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Bank  Angle,  degrees  c:  Angle  of  Attack,  degrees 
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Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 
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Angle  of  Attack  Profile,  Minimum  Downrange  Trajectories 

Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 
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Figure  57  Bank  Angle  Profile,  Minimum  Downrange  Trajectories 
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h 


Pseudo  Controls 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  58  Angle  of  Attaek  Rate  Profile,  Minimum  Downrange  Trajeetories 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  59  Bank  Angle  Profile,  Minimum  Downrange  Trajectories 


79 


L  Trim  Deficiency 


T rajectory  Trim  Deficiency 


Figure  60  Trajectory  Trim  Deficiencies,  Minimum  Downrange  Trajectories 


The  trim  deficiency  path  constraint  for  the  failed  vehicle  is  not  active  until 
passing  Mach  2  at  which  point  the  vehicle  follows  the  path  constraint  until  subsonic.  The 
initial  Mach-alpha  profiles  for  both  vehicles  are  nearly  identical  until  Mach  6  at  which 
point  they  diverge.  The  vehicles’  Mach-altitude  profiles  begin  to  slightly  diverge  until 
Mach  4  at  which  point  the  nominal  vehicle  conducts  a  climb.  While  it  appears  to  be 
some  encroachment  upon  the  trim  deficiency  path  constraint  for  the  failed  trajectory 
between  Mach  2  and  Mach  1,  all  trim  deficiency  points  corresponding  to  the  trajectory’s 
nodal  points  are  within  limits.  The  possible  encroachment  illustrates  the  need  to  have 
sufficient  nodes  not  only  for  endpoint  accuracy  but  also  to  adequately  characterize  the 
trajectory  as  related  to  trim  deficiency.  The  solution  can  avoid  encroachment  by  placing 
knots  at  the  critical  points  where  there  is  a  large  change  in  trim  deficiency  for  small 
changes  in  either  Mach  number  or  angle  of  attack. 
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j.  Optimality  Analysis 

The  trajectories  generated  by  DIDO  were  “near-optimal”.  An  examination 
of  the  Hamiltonians  for  both  trajectories  shows  larger  Hamiltonians  than  the  expected 
zero,  however,  the  average  Hamiltonian  for  both  trajectories  is  nearly  zero.  Both 
vehicles’  flight  path  angles  reach  the  limiting  value  of  ninety  degrees.  The  duals  for  both 
FPAs  meet  the  necessary  conditions. 


Hamiltonian  for  Nominal  Minimum  Downrange  Trajectory,  80  nodes 


Figure  61  Hamiltonian  for  Nominal  Minimum  Downrange  Trajectory 
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flight  path  angleDual 


States 


Figure  62  Constrained  Flight  Path  Angle  and  Associated  Dual,  Nominal  Minimum 

Downrange  Trajectory 


controls  dual 


controls 


Figure  63  Constrained  Angle  of  Attack  Rate  and  Associated  Dual,  Nominal 

Minimum  Downrange  Trajectory 
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Scaled  flight  path  angle  Dual  ^  Hamiltonian 


Hamiltonian  for  Stuck  Flap  Minimum  Downrange  Trajectory,  80  nodes 


64  Hamiltonian  for  Stuck  Flap,  Minimum  Downrange  Trajectory 


flight  path  angleDual 


States 


Figure  65  Constrained  Flight  Path  Angle  and  Assoeiated  Dual,  Stuck  Flap  Minimum 

Downrange  Trajeetory 
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controls  dual 


controls 


Figure  66  Constrained  Angle  of  Attack  Rate  and  Associated  Dual,  Stuck  Flap 

Minimum  Downrange  Trajectory 

k.  Feasibility 

Both  nominal  and  stuck  flap  trajectories  propagate  with  acceptable  error 
with  80-node  solutions.  Propagation  errors  are  identical  to  those  for  the  failed,  maximum 
downrange  trajectory. 
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DIDO  and  Propagated  States  for  Stuck  Flap  Trajectory,  80  nodes 


time,  TU 


Figure  67  DIDO  and  Propagated  Controls  for  Stuek  Flap,  Minimum  Downrange 
Trajectory  (DIDO  controls  are  denoted  by  discrete  points) 

5,  Minimum  /Maximum  Crossrange 

The  minimum  and  maximum  crossrange  profiles  show  symmetry  and  as  such 
only  the  maximum  crossrange  will  be  shown  as  representative  trajectories. 

a,  Mach-Altitude  Profile 

The  Mach-altitude  profile  is  similar  to  those  shown  for  the  maximum 
downrange  case.  The  vehicle  trades  airspeed  for  altitude  as  it  seeks  to  maximize  the  cross 
range. 
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Mach  -  Altitude  Profile  for  80-Node  T rajectory 


Figure  68  Downrange  Distance  Profiles  for  Maximum  Crossrange  Trajectories 
b.  Trajectory  in  Cartesian  Coordinates 

An  interesting  result  of  Figure  69  is  that  the  nominal  vehicle  reaches  the 
maximum  downrange  in  the  first  fifth  of  its  flight  time.  This  corresponds  to  an  azimuth 
of  90  degrees  at  the  300  second  point.  This  failed  vehicle  initiates  its  turn  at  the  same 
time  as  the  nominal  vehicle  but  is  able  to  reach  90  degrees  azimuth  (with  some  trouble) 
much  earlier  because  of  the  energy  bleed  caused  by  the  additional  drag  of  the  stuck  flaps. 
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Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  69  Downrange  Distance  Profile  for  Maximum  Crossrange  Trajectories 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  70  Crossrange  Distance  Profile  for  Maximum  Crossrange  Trajectories 
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Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  71  Altitude  Profile  for  Maximum  Crossrange  Trajeetories 

c.  Velocity  Profile 

Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  72  Velocity  Profile  for  Maximum  Crossrange  Trajectories 
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d.  Flight  Path  and  Azimuth  Angles 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  73  Flight  Path  Angle  Profile  for  Maximum  Crossrange  Trajeetories 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  74  Azimuth  Angle  Profile  for  Maximum  Crossrange  Trajeetories 
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e. 


Vertical  Speed  Profiles 


Nominal  Trajectory 


Trajectory  with  Flaps  Stuck  at  26  degrees 


Figure  75  Vertical  Speed  Profiles  for  Maximum  Crossrange  Trajectories 
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f.  Normal  Force  Profiles 


Normal  Force  for  Nominal  Max  Crossrange,  80-nodes 


Figure  76  Normal  Force  ( “  1 )  foi"  Nominal  Maximum  Crossrange  Trajectory 


Normal  Force  for  Failed  Max  Crossrange,  80-nodes 


Figure  77  Normal  Force( n^-l)  for  Stuck  Flap,  Maximum  Crossrange  Trajectory 
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Angle  of  Attack,  degrees 


g- 


Controls 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  78  Angle  of  Attack  Profile  for  Maximum  Crossrange  Trajectories 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  79  Bank  Angle  Profile  for  Maximum  Crossrange  Trajectories 
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h.  Pseudo  Controls 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  80  Angle  of  Attack  Rate  Profile  for  Maximum  Crossrange  Trajectories 


Comparison  of  Nominal  and  Stuck  Flap  Trajectories,  80-nodes 


Figure  81  Bank  Angle  Rate  Profiles  for  Maximum  Crossrange  Trajectories 
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T rajectory  Trim  Deficiency 


Angle  of  Attack,  Degrees 


Figure  82  Trim  Deficiencies  for  Maximum  Crossrange  Trajectories 


94 


Hamiltonian 


7- 


Optimality  Analysis 


Figure  83  Hamiltonian  for  Stuck  Flap,  Maximum  Crossrange  Trajectory 


Hamiltonian  for  80  node  Nominal  T rajectory 


Figure  84  Hamiltonian  for  Nominal  Maximum  Crossrange  Trajectory 
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Figure  85  Bank  Angle  Rate  and  Dual  for  Stuek  Flap  Max  Crossrange  Trajeetory 
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k.  Feasibility 


DIDO  and  Propagated  States  for  Norn  Trajectory,  80  Nodes 
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Figure  86  DIDO  and  Propagated  Controls  for  Stuek  Flap  Max  Crossrange  Trajeetory 
(DIDO  eontrols  are  denoted  by  discrete  points) 


The  error  between  the  DIDO  states  and  propagated  states  for  the  vehicle 
with  stuck  flaps  is  shown  below. 


error  = 


37/t 


-85/t  -1ft  -7  deg  9  deg  -0.01^^^  -0.01 

■>  ■>  /sec  ^  ^  /sec  /  sec 


C.  TRAJECTORY  TO  A  SPECIFIED  FIELD 

1,  Endpoint  Conditions  and  Cost  Function 

The  endpoint  set  for  a  trajectory  to  a  specified  field  within  the  footprint  contains 
the  field  Cartesian  coordinates  as  well  as  the  azimuth  of  the  landing  runway.  Both 
nominal  and  failed  vehicles  were  able  to  land  at  specified  fields  within  their  respective 
footprints  however  the  choice  of  runway  azimuths  was  limited  to  an  approximate  range 
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of  +/-  100  degrees.  For  example,  landing  at  a  runway  within  the  failed  footprint  with  a 
runway  heading  of  -90  degrees  was  aeeomplished  relatively  rapidly  with  a  resultant 
“loeally  optimal  solution  found”.  The  same  field  with  a  runway  heading  of  -100  degrees 
gave  a  “near  optimal”  solution  while  a  field  with  a  runway  heading  of  -120  degrees  gave 
infeasible  or  suboptimal  results  for  the  same  80-node  run.  It  should  be  noted  that  while 
the  infeasibility  was  for  an  80-node  trajeetory,  larger  node  trajeetories  have  been  shown 
to  provide  near  optimal  or  loeally  optimal  solutions  when  bootstrapped  from  lower-node 
infeasible  trajeetories.  There  is  the  additional  eomputation  time  penalty  that  eomes  along 
with  eomputing  trajeetories  with  larger  nodes. 

The  eost  funetion  for  these  trajeetories  to  speeified  field  was  a  Lagrangian  eost 
funetion  that  minimized  the  angles  of  attaek  and  bank.  The  rationale  behind  a  eost 
funetion  sueh  as  this  is  that  depending  on  the  eontrol  system  failure,  exeessive  movement 
of  eontrol  surfaee  hydraulie  aetuators  may  deplete  the  hydraulie  working  fluid  of  that 
partieular  loop  resulting  in  easeading  failures. 

The  full  problem  formulation  is  shown  below: 

min„  T(x(t),M =  j 

h) 

subjeet  to  x  =  f  (x,  u,  r) 
hi  - 

Where 

=  yf  Zf 

=  [110  miles  -1  mile  500  feet  -500 ±1000  ft/ min  Mach  0.15±2%  -100°  0±10°]'' 
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All  other  initial  conditions  and  path,  state  and  control  vectors  and  limits  are 
identical  to  the  previous  formulations. 

2,  Profiles  for  a  Failure-Mode  Trajectory  to  a  Specified  Field 
a.  Mach-Altitude  Profile 


Mach  -  Altitude  Profile  for  80-Node  Trajectory 


Figure  87  Mach-altitude  Profile  for  Stuck  Flap  Trajectory  to  Specified  Field 
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Downrange  Distance,  miles 


b.  Trajectory  in  Cartesian  Coordinates 

80-node  T rajectory  to  Specific  Field 


Figure  88  Downrange  Profile  for  Stuck  Flap  Trajectory  to  Specified  Field  (Flaps 

Stuck  at  26  degrees) 
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80-node  T rajectory  to  Specific  Fieid 


Figure  89  Crossrange  Profile  for  Stuck  Flap  Trajectory  to  Specified  Field  (Flaps 

Stuck  at  26  degrees) 


time,  seconds 


Figure  90  Altitude  Profile  for  Stuck  Flap  Trajectory  to  Specified  Field 
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3-Dimension  Piot  of  Entry  T rajectory 


Figure  91  Three  Dimensional  Trajeetory  Depletion  (not  to  seale)  for  Failed  Vehiele 

(Flaps  Stuek  at  26  degrees) 
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Altitude,  miles 


Figure  92  Quiver  Plot  (not  to  scale)  for  Failed  Vehicle  Trajectory  to  Specified  Field 

(Flaps  Stuck  at  26  degrees) 
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Ground  Track  Plot  of  Trajectory 


Figure  93  Ground  Track  Plot  (not  to  scale)  of  Failed  Vehicle  Trajectory  to  Specified 

Field  (Flaps  Stuck  at  26  degrees) 
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Velocity,  mph 


C, 


Velocity  Profile 


80-node  T rajectory  to  Specific  Field 


Figure  94  Velocity  Profile  for  Failed  Vehicle  Trajectory  to  Specified  Field 
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d.  Flight  Path  and  Azimuth  Angle 


80-node  T rajectory  to  Specific  Fieid 


Figure  95  Flight  Path  Angle  Profile  for  Failed  Vehiele  Trajeetory  to  a  Speeified 


Field 

80-node  T rajectory  to  Specific  Fieid 


Figure  96  Azimuth  Angle  Profile  for  Failed  Vehicle  Trajectory  to  a  Specified  Field 
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e. 


Body  Frame  Normal  Force 


stuck  Body  Flap,  80  node  T rajectory 


Figure  97  Normal  Force(  n^-l)  and  Angle  of  Attack  Profile  for  Failed  Vehicle 

Trajectory  to  Specified  Field 
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stuck  Body  Flap,  80  node  T rajectory 


Figure  98  Normal  Force(  n^-l)  and  Angle  of  Bank  Profile  for  Failed  Vehicle 

Trajectory  to  Specified  Field 
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f.  Vertical  Speed  Profile 


Vertical  Speed  For  26  deg  Stuck  Flap  Trajectory  to  Specific  Field 


distance,  miles 

Figure  99  Vertical  Speed  Profile  for  Failed  Vehicle  Trajectory  to  Specified  Field 
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Vertical  Speed  For  26  deg  Stuck  Flap  Trajectory  to  Specific  Field 


Figure  100  Vertical  Speed  Profile  (vs.  time)  for  Failed  Vehicle  Trajectory  to  Specified 

Field 
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80-node  T rajectory  to  Specific  Fieid 


Figure  101  Angle  of  Attaek  Profile  for  Failed  Vehiele  Trajeetory  to  Speeified  Field 


80-node  T rajectory  to  Specific  Fieid 


Figure  102  Angle  of  Bank  Profile  for  Failed  Vehicle  Trajectory  to  Specified  Field 
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k 


Pseudo  Controls 


80-node  Trajectory  to  Specific  Fieid 


Figure  103  Angle  of  Attack  Rate  Profile  for  Failed  Vehicle  Trajectory  to  Specified 

Field 
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Bank  Angle  Rate,  degrees/second 


The  bank  angle  rate  profile  is  nearly  bang-bang  for  this  trajectory. 
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Figure  104  Bank  Angle  Rate  for  Failed  Vehicle  Trajectory  to  Specified  Field 
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u 


Trim  Deficiency 


T rajectory  Trim  Deficiency 


Figure  105  Trim  Deficiency  Profile  for  Failed  Vehicle  Trajectory  to  Specified  Field 


The  trim  deficiency  profile  illustrates  several  problems  associated  with 
fitting  curves  to  the  trim  deficiency  data.  The  contour  plot  of  the  fitted  trim  deficiency 
curve  is  relatively  well  behaved  except  for  the  low  Mach  and  angle  of  attack 
combinations,  which  is  exactly  where  the  trajectory  wants  to  go.  For  this  case,  the 
process  of  fitting  the  trim  deficiency  data  has  artificially  reduced  the  Mach-alpha 
envelope  available  to  the  trajectory  with  a  corresponding  decrease  in  the  performance  of 
the  vehicle.  There  must  be  sufficient  nodes  or  properly  placed  knots  as  previously 
discussed  to  cover  regions  in  the  trajectory  where  there  is  rapidly  changing  trim 
deficiencies.  While  the  solutions  points  do  not  violate  the  trim  deficiency  path  constraint 
as  shown  below,  the  trajectory  between  the  points  may  violate  the  constraint.  This 
violation  can  be  easily  removed  with  proper  problem  setup. 
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Columns  1  through  15 


0.0005  0.0004  0.0004 

0.0003 

0.0002 

0.0001 

0.0000 

0.0001 

0.0001 

0.0002 

0.0004 

0.0005 

0.0006 

0.0003  ■ 

0.0008 

Columns  16  through  30 
-0.0008  0.0002  0.0005 

0.0004 

0.0001 

-0.0002 

0.0001 

0.0014 

0.0025 

0.0030 

0.0030 

0.0030 

0.0024 

0.0030 

0.0030 

Columns  3 1  through  45 
0.0021  0.0022  0.0030 

0.0030 

0.0030 

0.0030 

-0.0020 

-0.0033 

-0.0013 

0.0011 

0.0030 

0.0030 

0.0011 

-0.0025 

-0.0081 

Columns  46  through  60 
-0.0145  -0.0182  -0.0201 

-0.0278 

-0.0412 

-0.0591 

-0.0806 

-0.1083 

-0.1281 

-0.1283 

-0.1106 

-0.0766 

-0.0514 

-0.0352 

:  -0.0412 

Columns  61  through  75 
-0.0585  -0.0721  -0.0858 

-0.0912 

-0.0861 

-0.0600 

0.0002 

0.0001 

-0.0000 

-0.0001 

-0.0002 

-0.0001 

-0.0002 

-0.0001 

-0.0001 

Columns  76  through  80 
-0.0000  0.0000  0.0001  0.0001  0.0001 

Table  2.  Numerical  Trim  Deficiency  Constraint  Results  for  80  nodal  Points 


j.  Optimality  Analysis 


stuck  Body  Flap,  80  node  T rajectory 


Figure  106  Flamiltonian  for  Failed  Vehicle  Trajectory  to  Specified  Field  (Flaps  Stuck 

at  26  degrees) 
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Scaled  flight  path  angle  Dual 


Figure  107 


flight  path  angleDual 


FPA  and  Associated  Dual  for  Failed  Vehicle  Trajectory  to  Specified  Field 
(Flaps  Stuck  at  26  degrees) 
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Scaled  bank  angle  rate 


controls  dual 


Figure  108  Bank  Angle  Rate  and  Associated  Dual  for  Failed  Vehicle  Trajectory  to 
Specified  Field  (Flaps  Stuck  at  26  degrees) 
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Scaled  aoa  rate 


controls  dual 


Figure  109  Angle  of  Attack  Rate  and  Associated  Dual  for  Failed  Vehicle  Trajectory 
to  Specified  Field  (Flaps  Stuck  at  26  degrees) 
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k.  Feasibility  Analysis 


DIDO  and  Propagated  States  for  Stuck  Flap  Trajectory,  80  Nodes 


Figure  110  Scaled  DIDO  and  Propagated  States,  Failed  Vehicle  Trajectory  to 
Specified  Field  (DIDO  controls  are  denoted  by  discrete  points) 


The  error  between  the  DIDO  states  and  the  R-K  propagated  states  using 
the  DIDO  controls  is  shown  in  Figure  XX.  The  altitude  error  is  still  high  and  either  the 
trajectory  should  be  generated  using  more  nodes  or  the  endpoint  conditions  relaxed 

Error  =  [108ft  1247ft  -460ft  -01.7ft/sec  -O.Sdeg  4.5deg  Odeg/sec  -O.ldeg/sec] 
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D,  DIDO  PERFORMANCE 

While  there  are  several  performanee  eriteria  for  the  optimal  trajeetory  solutions, 
several  eriteria  dominate; 

i.  Aecuraey  of  the  solution.  In  other  words,  how  closely  do  the  states  at  the 
endpoint  set  generated  by  using  the  propagated  solution  controls  match  the 
internally  generated  states  used  by  the  optimization  algorithm. 

ii.  Optimality.  How  optimal  is  the  solution?  Examination  of  the  necessary 
conditions  can  help  answer  this  question  but  can  be  inadequate  when 
attempting  to  qualify  a  “near  optimal”  solution. 

iii.  Time  to  find  a  solution.  Can  the  solution  be  found  in  “real  time”  and  if 
not,  how  close  to  real  time  can  a  solution  be  found? 

For  the  psuedospectral  method  used  in  this  thesis,  all  three  criteria  are  a  function 
of  the  number  of  nodes  used  in  the  solution.  Given  a  feasible  solution,  the  error  between 
propagated  and  DIDO  states  can  be  reduced  by  increasing  the  number  of  nodes.  At  the 
same  time,  increasing  the  number  of  nodes  in  the  solution  comes  at  significant  expense  in 
computation  time.  If  a  locally  optimal  solution  is  declared,  the  end  states  will  not 
generally  change;  however,  if  a  “near  optimal”  solution  is  declared,  then  the  end  states 
become  a  function  of  the  number  of  nodes  in  the  solution.  Recall  that  the  generation  of  a 
“near  optimal”  solution  is  a  consequence  of  many  things  including  insufficient  number  of 
nodes  and  the  use  of  nonsmooth  data. 

The  following  analysis  provides  some  insight  into  the  DIDO  solution  performance 
for  a  nominal  trajectory  with  the  same  initial  and  endpoint  conditions  as  shown  in  the 
results.  The  graphs  represent  data  from  a  nominal  maximum  downrange  trajectory.  The 
relationship  between  propagated  errors  and  number  of  nodes  in  the  solution  is  shown  in 
Figure  III  and  Figure  1 12  Notice  that  the  error  for  low  node  solutions  is  unacceptable 
while  the  error  for  solutions  greater  than  60  nodes  results  is  acceptable. 
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error,  miles 


Error  Between  DIDO  and  Propagated  States 


Figure  111  Downrange  and  Crossrange  Error,  Nominal  Trajectory 


Error  Between  DIDO  and  Propagated  States 


Figure  1 12  Altitude  Error  vs.  Nodes 
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DIDO  execution  times  for  P4  2.0  Ghz  Processor 


Figure  113  Time  to  Reaeh  a  Solution,  Nominal  Trajeetory  for  Pentium  4,  2.0  Ghz 

processor  using  Windows  2000. 

While  larger  nodes  give  better  propagation  performance,  the  cost  is  significant. 
The  execution  time  for  an  80-node  solution  is  nearly  double  that  of  a  60-node  solution. 
These  times  were  for  unoptimized  code  running  on  a  networked  Pentium  4  processor 
running  Matlab  6.0  in  a  full  Windows  2000  environment.  Considerable  performance 
increases  resulting  from  speed-optimized  code  and  a  faster  processor  using  a  dedicated 
(non-multitasking)  computing  environment  can  be  expected.  Execution  times  of  70 
seconds  for  20-  node  nominal  trajectories  have  been  routinely  observed  on  a  P4  3.06  Ghz 
laptop  with  minimum  cpu  overhead. 

Figure  114  and  Figure  115  illustrate  the  convergence  of  the  “near  optimaf’ 
solution.  Notice  that  by  60-nodes  the  final  DIDO  cross  and  downrange  state  variables  are 
nearly  identical  while  altitude  convergence  is  reached  immediately.  The  intent  of  these 
graphs  is  to  show  the  pitfall  of  choosing  a  small  number  nodes  to  generate  faster 
solutions.  The  difference  between  the  maximum  downrange  of  a  20  node  solution  and  a 
60  node  solution  is  nearly  ten  miles. 


122 


Convergence  of  Near  Optimal  Solutions 


Figure  1 14  Convergence  of  Down  and  Cross  Ranges  for  Near  Optimal  Solutions 
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Altitude  Convergence  for  Nominal  Near  Optimal  Trajectory 
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V.  FUTURE  WORK 


A.  CONVERSION  TO  SI  UNITS 

This  thesis  was  conducted  using  the  English  unit  system.  The  units  for  velocity, 
heating  rate  and  pressure  do  not  lend  themselves  to  convenient  balancing  or  scaling. 
Ross  [ref  21]  has  shown  that  the  use  of  English  units  in  the  problem  formulation  can 
result  in  numerical  difficulties  leading  to  infeasible  solutions.  The  metric  system  is  a 
more  intuitive  system  with  all  units  a  subset  of  the  primary  distance,  mass  and  time  units 
of  meters,  kilograms  and  seconds.  Problems  formulated  using  the  inks  system  lead  to 
more  easily  balanced  and  scaled  units  with  better  corresponding  numerical  results. 


B,  DATA  SMOOTHING 

The  DIDO  software  package  performance  is  best  with  smoothed  data  that  avoids 
singularities  in  the  gradients  of  the  state  variables.  The  presence  of  non  smooth  data 
results  in  extremely  long  calculation  times,  “near  optimal”  vice  “locally  optimal” 
solutions  along  with  the  presence  of  “sweet  spots”  where  a  small  range  of  nodes  results  in 
feasible  solutions.  One  way  to  artificially  smooth  the  data  is  to  increase  the  number  of 
nodes  in  the  solution;  however,  the  computation  penalty  can  be  substantial.  The  best 
approach  may  be  to  use  a  more  finely  meshed  lookup  Cd  and  Cl  table  in  the  dynamics 
file.  These  tables  are  generated  offline  and  have  minor  impact  on  the  lookup  routine;  the 
extra  lookup  time  is  easily  offset  by  the  quicker  convergence  time. 


C.  TRIM  DEFICIENCY  MAPPING 

The  trim  deficiency  function  used  in  this  thesis  was  a  polynomial  fitted  to  the 
tabular  trim  deficiency  data.  This  was  done  to  minimize  execution  time  of  the  solution. 
The  curve  fit  to  the  data  was  not  perfect  and  a  better  approach  may  be  to  use  a  finely 
meshed  lookup  table  as  discussed  in  the  previous  section.  This  would  help  automate  the 
trim  deficiency  process  with  some  penalty  in  execution  time. 
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D,  CLOSED  LOOP  IMPLEMENTATION 

While  80  node  solutions  ean  take  up  to  50  minutes  to  eonverge  on  a  Pentium-4 
2.0  ghz  proeessor,  an  optimized  eode  using  40  nodes  ean  aehieve  eonvergenee  in  as  little 
as  1  minute  on  the  same  maehine.  While  40-node  solutions  do  not  generally  provide  the 
neeessary  aeeuraey  for  the  endpoints,  it  may  be  possible  to  reduee  outer  loop  exeeution 
time  using  lower  order  solutions  and  still  meet  the  aeeuraey  requirements  by  using 
feedbaek.  This  should  be  more  elosely  investigated  to  determine  the  feasibility  of  this 
approaeh. 


E.  OPTIMIZING  THE  DIDO  CODE  FOR  SPEED 

The  algorithms  used  in  this  thesis  were  eonstrueted  with  minimal  emphasis  on 
eomputation  speed.  There  are  many  options  available  to  inerease  the  speed  of  the 
algorithm.  The  first  and  most  signifieant  speed  inerease  ean  be  seen  by  removing  the 
eode  from  the  Windows  (multi-tasking)  environment  and  eompiling  for  use  on  a 
dedieated  real  time  proeessor.  Additional  improvements  in  the  MATLAB  eode  ean 
provide  further  speed  inereases;  MATLAB  7  shows  some  speed  inerease  over  MATLAB 
6.5.  As  CPU  speed  eontinues  to  inerease  and  the  operating  eode  is  optimized,  exeeution 
times  should  eontinue  to  deerease.  Additional  future  improvements  to  the  DIDO 
software  paekage  have  the  potential  of  aehieving  signifieant  deereases  in  solution  times. 


F.  FULL  6-DOF  MODEL  WITH  INNER  AND  OUTER  LOOP  COUPLING 

As  detailed  in  the  footprint  diseussion,  the  deeoupled  inner  and  outer  loops  along 
with  a  redueed  order  model  do  a  satisfaetory  job  of  eomputing  the  vehiele  performanee. 
The  true  performanee  of  the  vehiele;  however,  is  not  realized  beeause  the  two  loops  are 
deeoupled  as  illustrated  in  the  differenees  in  the  nominal  and  failed  vehiele  trajeetory 
footprints.  While  the  use  of  a  6-DOF  model  with  inner  and  outer  loop  may  appear  to  be 
the  more  diffieult  problem  to  solve,  the  eomplex  interplay  between  the  dynamieal 
equations  and  state  and  path  eonstraints  may  reduee  the  problem  to  one  more  easily 
solved.  The  potential  benefits  of  having  a  fully  integrated  inner  and  outer  loop  are  great 
enough  to  warrant  further  investigation  into  the  feasibility  of  this  approaeh. 
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